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CHAPTER  1:  SPACE  PROPULSION  SYSTEMS 


1.  J  Rocket  Propulsion  Fundamentals 

A  rocket  is:  a  device  which  propels  a  vehicle  by  expelling  mass.  In 
space,  where  aerodynamic  forces  arc  negligible  and  where  solar  radiation  pressure 
forces  are  assumed  to  be  negligible,  the  equation  of  motion  of  a  rocket  powered 
vehicle  is  given  by  Eq.  (1. 1). 

m  v  ~  me  +  mg 

In  this  equation  m  is  the  mass  of  the  vehicle,  r  is  the  position  vector,  c  is  the 
exhaust  velocity  vector  of  the  rocket,  g  is  the  gravitational  acceleration,  and  m 
is  the  rate  at  which  mass  is  expelled  from  the  rocket.  The  rocket  exhaust  is 
assumed  to  be  uniform  in  direction  and  magnitude.  This  equation  gives  the  proper 
form  of  Newton's  law  for  a  variable  mass  body  whose  loss  of  mass  is  solely 
through  a  uniform  rocket  exhaust  stream.  The  thrust  of  the  rocket  is  defined 
by  Eq.  (1.2), 

f  =  m  c , 

and  since  the  rate  at  which  the  vehicle  loses  mass  is  negative,  the  thrust  will  be 
directed  opposite  to  the  exhaust  velocity  vector  of  the  rocket. 

For  electrically  propelled  vehicles,  a  quantity  of  some  importance  is 
the  power  in  the  exhaust  stream  since  this  power  must  be  supplied  from  a  power 


source.  The  exhaust  stream  power  is  given  by  Eq.  (1.3)  for  a  uniform  exhaust 
stream. 


(1.3) 


The  rate  of  mass  flow  is  proportional  to  the  first  power  of  thrust  for  a  rocket 
with  constant  exhaust  velocity  and  is  proportional  to  the  second  power  of  thrust 
for  a  rocket  with  constant  power.  It  is  inversely  proportional  to  either  exhaust 
velocity  or  power,  hence  it  is  desirable  to  have  high  exhaust  velocity  or  high 
power  in  rockets.  Instead  of  exhaust  velocity,  rocket  engineers  normally  speak 
of  a  quantity  called  specific  impulse  which  is  most  conveniently  defined  as  the 
exhaust  speed  divided  by  the  standard  acceleration  of  gravity,  Kq.  (1.5). 


Specific  impulse  is  normally  measured  in  units  of  seconds. 

The  variation  of  thrust  with  mass  flow  rate  for  the  two  types  of  rockets 
is  shown  in  Figure  1. 1.  It  should  be  noted  that  the  curves  cross  at  two  points. 
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.1.1  Variation  of  Thrust  With 
Mass  Flow  Rate 


Constant  Exhaust  Velocity 


at  the  origin  and  at  some  finite  thrust.  Many  rockets  cannot  be  operated  at 

/ 

various  thrust  levels  but  can  be  switched  on  and  off.  Such  rockets  may  be  con¬ 
sidered  as  having  either  the  same  power  or  the  same  exhaust  velocity  at  their 
operating  point  and  when  they  are  switched  off. 

1.2  The  Spectrum  of  Space  Propulsion  Systems 

Any  rocket  engine  may  be  characterized  by  two  important  quantities 
which  affect  vehicle  performance:  the  ratio  of  thi'ust  to  engine  weight  for  the 
rocket  and  the  specific  impulse  of  the  rocket.  The  ratio  of  thrust  to  engine 
weight  is  a  measure  of  the  mass  of  the  engine  required  for  a  given  thrust  level, 
while  the  specific  impulse  gives  a  measure  of  the  fuel  required  for  various  mis¬ 
sions.  The  possible  variations  in  these  quantities  for  space  propulsion  systems 
are  surprisingly  wide,  falling  basically  into  two  different  areas.  A  plot  of  the 
ratio  of  the  thrust  to  engine  weight  (at  standard  gravity)  versus  specific  impulse 
is  shown  in  Figure  1.2  for  some  of  the  more  conventional  types  of  rockets.  All 
the  way  up  in  the  left-hand  corner  of  the  figure,  there  is  a  block  representing 
chemical  rockets.  Chemical  rockets  liavc  large  ratios  of  thrust  to  engine  weight 
The  engines  weigh  very  little  for  the  thrust  they  produce,  but  they  also  have  rela 
tively  low  specific  impulses.  There  are  three  major  types  of  chemical  rockets: 
those  using  liquid  propellants  which  are  stored  in  a  tank  and  burned  in  a  rocket 
engine,  those  which  use  solid  propellants  where  the  propellant  lank  is  also  the 
thrust  chamber  of  the  engine,  and  those  using  both  liquids  and  solids  where  the 


( 
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Fig.  1.2  Space  Propulsion  Engines 


/■ 


solid  propellant  is  again  stored  in  the  thrust  chamber  while  the  liquid  propellant 
is  stored  in  a  separate  tank.  Liquid  propellants  tend  to  offer  higher  performance 
than  solid  propellants,  having  both  higher  specific  impulses  and  generally  lighter 
overall  system  weight  (including  both  tanks  and  engine).  Solid  propellants  do 
have  advantages  such  as  storability,  density,  and  economy  which  make  them 
desirable  for  some  missions.  Hybrid  rockets  using  both  solids  and  liquids  have 
not  been  developed  very  extensively  yet  but  do  give  promise  of  fairly  high  specific 
impulses  and  may  be  promising  for  some  missions.  Only  liquid  rockets  have  been 
shown  in  the  figure  as  it  does  not  include  tankage  weights. 

Liquid  rockets  used  today  with  hydrogen-oxygen  propellants  arc  pro¬ 
ducing  specific  impulses  of  approximately  430  seconds  and  some  more  advanced 
systems  have  been  proposed  which  may  produce  specific  impulses  on  the  order 
of  500  seconds.  However,  it  does  not  seem  likely  that  chemical  rockets  may 
produce  performance  much  in  excess  of  this  if  stable  chemical  species  are  used. 

Liquid  and  hybrid  chemical  rockets  may  be  designed  so  that  they  can 

''  / 

be  throttled.  In  this  case,  the  engine  can  have  any  thrust  varying  from  some 
maximum  thrust  down  to  some  fairly  low  value.  This  throttling  takes  place  essen¬ 
tially  at  constant  specific  impulse. 

Because  chemical  rockets  have  relatively  low  specific  impulse,  a 
number  of  higher  specific  impulse  devices  have  been  proposed  and  some  are  under 
development.  None  of  these,  however,  has  yet  reached  the  application  stage  for 
primary  propulsion.  One  such  scheme  is  the  solid-core  nuclear  rocket  shown 
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by  a  block  next  to  llic  chemical  rocket.  These  rockets  utilize  hydrogen  exhausted 
through  a  very  high  temperature  nuclear  reactor  operating  near  the  material 
limits  of  graphite  or  tungsten  and  can  produce  specific  impulses  on  the  order  of 
800  seconds,  about  twice  that  obtainable  with  chemical  rockets.  This  large  im¬ 
provement  in  specific  impulse  is  obtained  at  the  expense  of  some  increase  in 
both  engine  and  tankage  weights  for  this  type  of  rocket,  but  they  arc  advantageous 
for  many  missions  and  there  is  a  current  project  to  develop  a  solid-core  nuclear 
rocket.  Like  chemical  rockets,  nuclear  rockets  maybe  designed  to  be  throttled 
and  may  do  this  at  essentially  constant  specific  impulse. 

Both  chemical  rockets  and  solid-core  rockets  are  considered  to  be  high- 
thrust  devices  that  can  be  used  for  taking  off  from  the  surface  of  planets.  They 
can  provide  high  enough  accelerations  so  that  the  rocket  vehicle  performance  may 
often  be  approximated  by  replacing  the  finite  burning  time  of  these  rockets  by  an 
impulse  which  requires  zero  time.  Their  primary  limitation  is  exhaust  velocity. 
Because  of  their  throttling  characteristics,  they  will  be  referred  to  as  constant- 
exhaust-vclocity  rockets  throughout  this  book. 

Running  along  the  bottom  of  Figure  1.2  is  a  class  of  electric  propulsion 

systems  which  operate  at  much  higher  specific  impulses  than  the  high-thrust 

/ 

chemical  and  nuclear  rocket  systems,  but  also  have  orders  of  magnitude  larger 
engine  weights  for  a  given  thrust.  These  devices  can  only  be  used  in  oi'bit. 

Their  thrust  is  very  much  smaller  than  local  gravity  for  close  planetary  orbits 
and  their  whole  mode  of  operation  in  a  space  mission  is  quite  different.  In  a 


s 
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chemical  or  nuclear  rocket,  the  powered  flight  time  is  on  the  order  of  minutes 
or  hours  whereas  for  an  electric  propulsion  system  the  powered  flight  time  is 
on  the  order  of  months  or  years  and  generally  constitutes  the  majority  of  the 
total  mission  time.  An  electric  propulsion  system  is  a  combination  of  a  power 
source  and  suitable  power  conditioning  equipment  with  a  thrust  producing  device 
or  thrustor.  The  same  types  of  thrustors  might  be  used  with  different  types  of 
power  sources  and  vice  versa.  A  great  variety  of  thrustors  have  been  proposed 
and  some  of  these  have  been  developed  to  an  operational  state.  The  only  currently 
operational  thrustor  types  arc  resistojets  which  have  relatively  modest  specific 
impulses  of  up  to  800  seconds.  Another  highly  developed  type  of  Ihrustor  is  the 
ion  rocket  which  operates  efficiently  at  specific  impulses  above  3000  seconds.  A 
number  of  other  devices  such  as  nrejets,  MPD  arcs  and  various  plasma  devices 
have  been  proposed  and  some  of  these  are  under  development.  Most  of  these 
engines  can  be  operated  only  at  or  near  a  design  point  and  cannot  vary  thrust 
effectively.  It  would  be  desirable  to  develop  thrustors  which  could  operate  effi¬ 
ciently  over  a  wide  range  of  specific  impulses  at  essentially  constant  power. 

This  can  be  approximated  to  some  extent  by  carrying  several  thrustors  on  a 

I 

given  power  supply  and  switching  from  one  to  the  other. 

It  will  be  noted  that  the  thrust-to-weight  ratio  versus  specific  impulse 
for  these  electric  propulsion  devices  is  inversely  proportional  to  the  specific 
impulse  and  this  is  a  simple  consequence  of  Eq.  (1.3)  of  the  last  section.  The 
weight  of  an  electric  propulsion  system  is  primarily  in  the  power  supply  and  as 
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the  specific  impulse  is  raised,  the  thrust  must  be  decreased  in  order  to  keep  the 
power  constant.  Most  early  studies  of  electric  propulsion  considered  nuclear- 
electric.  power  supplies.  These  are  devices  which  take  the  heat  from  a  nuclear 
reactor  and  convert  it  in  some  fashion  or  another  to  electricity.  At  the  moment 
there  is  r.o  active  program  in  this  country  to  develop  a  nuclear-electric  power 
supply  which  would  be  suitable  for  electric  propulsion  and  attention  has  switched 
to  solar-electric  propulsion  systems.  There  are  many  possible  ways  of  converting 
solar  energy  into  electricity,  but  the  most  widely  used  and  the  most  successful 
c-f  these  has  been  the  photovoltaic  cell.  Unlike  the  nuclear-electric  systems,  the 
power  of  a  solar  cell  system  Is  a  function  of  distance  from  the  sun.  In  addition, 
the  solar  cell  array  must  be  in  full  sunlight  and  properly  oriented  towards  the  sun 
in  order  to  develop  full  power.  The  power  developed  is  not  an  inverse  square 
function  of  distance  from  the  sun  because  as  the  temperature  of  the  cells  changes, 
their  efficiency  varies.  A  typical  curve  for  power  as  a  function  of  distance  from 
the  sun  is  given  in  Figure  1.3. 

A  great  variety  of  other  propulsion  devices  has  been  proposed.  Among 
these  are  solar  heated  rockets,  radioisotope  rockets,  liquid  core  nuclear  rockets, 
gaseous  core  nuclear  rockets,  nuclear  pulse  rockets  and  solar  sails.  With  the 
exception  of  the  last  device,  these  are  all  rockets  and  may  be  treated  by  the  methods 
to  be  developed  in  this  book.  The  solar  sail  is  a  somewhat  academic  system  which 
is  easily  treated  by  the  general  principles  to  be  illustrated  in  this  monograph. 
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Distance  to  Sun  (  an) 
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1.  3  rayoff  Functions 


The  payoff,  the  thing  that  we  shall  bo  trying  to  minimize,  is  the  total 
fuel  used  and  is  given  by  integrating  Eq.  (1.4)  to  give  Eq.  (1.6), 

t1  t1 

m^  -  -  J  —  dt  -  J*  dt 

0  C  0  2P 

Two  forms  of  this  equation  are  given:  one  applicable  to  systems  with  constant 
exhaust  velocity  and  one  applicable  to  the  power-limited  electric  propulsion 
systems.  For  some  purposes  it  is  convenient  to  use  alternate  forms  for  these 
payoff  functions,  particularly  for  the  idealized  cases  where  the  engines  will  have 
variable  thrust.  If  we  have  a  constant  exhaust  velocity  system,  Eq.  (1.4)  may  be 
integrated  by  dividing  through  by  the  mass  to  form  Eq.  (1.  7). 

0  ^  f 

c  l  n  ~  =  T  -  dt  s  AV 
m1  6  m 

The  logarithm  of  the  mass  ratio  of  the  system,  a  monotonic  function  of  the  fuel 
used,  is  proportional  to  the  time  integral  of  the  thrust  acceleration  of  the  vehicle. 
This  quantity  is  spoken  of  as  characteristic  speed,  or  velocity  increment,  or  by  a 
great  variety  of  other  names  and  will  be  used  frequently.  It  is  possible  to  use 
this  payoff  function  instead  of  the  fuel  used  as  it  is  a  monotonic  function  of  the 
fuel  used.  The  characteristic  speed  may  be  interpreted  physically  as  the  change 
in  velocity  lint  would  be  produced  by  thrusting  in  a  single  direction  in  the  absence 
of  a  gravity  field.  Without  gravity,  the  change  in  speed  would  be  independent  of 
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the  acceleration  program.  Alternately,  the  characteristic  speed  may  be  looked 
upon  as  the  sum  oi  the  absolute  magnitude  of  the  velocity  increments  that  would 
be  produced  by  impulsive  thrusting  in  a  gravity  field. 

It  should  be  noted  that  the  characteristic  speed  can  be  made  much  larger 
than  tlic  exhaust  velocity.  It  the  final  mass  could  approach  zero,  the  character¬ 
istic  speed  would  approach  infinity. 

For  n  power-limited  system,  the  corresponding  quantity  is  found  by 
dividing  Eq.  (1.4)  through  by  the  square  of  the  mass  to  produce  Eq.  (1.8). 


p!.  rl  i  pH.  / 1  \2  dt  . , 


This  integral  lias  never  acquired  a  popular  name  and  is  generally  known  as  J. 
The  important  thing  here  is  that  another  monotonic  function  of  the  mass  of  the 
vehicle  is  proportional  to  the  square  of  the  thrust  acceleration  so  that  for  a 
power-limited  vehicle  it  is  desirable  to  minimize  the  lime  integral  of  the  square 
of  the  acceleration,  whereas  for  a  constant  exhaust  velocity  system,  it  is  desir¬ 
able  to  minimize  the  time  integral  of  the  acceleration  itself. 

It  was  mentioned  in  the  first  section  tliat  a  rocket  that  can  be  switched 
on  and  off  and  has  only  one  operating  point  can  be  regarded  as  either  a  constant 
power  or  a  constant  exhaust  velocity  rocket.  Either  type  of  payoff  may  be  used 
for  such  a  rocket.  In  fact,  in  this  case  there  is  the  following  simple  relation 
between  the  two  payoff  variables. 
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J 


<1.9) 


0 


m 


AV/c  ,  "I 

c  -l] 


1.4  Mass  Relations  for  Kxlinust-Volocity-I.imilccl  Systems 

The  total  mass  of  the  space  vehicle  is  assumed  to  consist  of  four  types 
of  muss:  payload,  powerplant,  structure  and  fuel  (or  propellant) 


0 

m  =  mL  +  mp  +  m  +  rn^ 


(1.  10, 


Tlte  powerplant  mass  and  part  of  the  structure  mass  is  assumed  to  be  proportional 
to  the  maximum  engine  thrust. 


mp  +  msp 


ki  f 

1  max 


<1.11) 


The  remainder  of  the  structure  mass  is  assumed  to  be  proportional  to  the  mass 
of  propellant. 


"sf  "  k2  mr 


(1.12) 


The  ratio  of  the  payload  mass  to  the  initial  mass  is  then  given  by  Eq.  (1. 13). 


mL  f  mF 

-T  “ 

m  m  m 


(1.13) 
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By  using  Eq.  (1.7)  of  the  last  section,  Eq.  (1.13)  may  be  put  into  tbc  desired  form, 


vehicle  is  obtained. 


given  value  for  exhaust  beam  power  and  powerplant  mass.  In  this  latter  case, 
and  only  in  this  latter  case,  it  is  possible  to  optimize  the  mass  distribution  in 
the  stage  independently  of  the  trajectory. 

For  the  idealized  power-limited  rocket,  having  fully  variable  thrust, 
the  value  of  J  in  Eq.  (1. 17)  will  depend  only  upon  the  particular  mission  and  will 
be  independent  of  the  vehicle  design.  If  the  relationship  between  powerplant  mass 
and  exhaust-beam  power  is  given  by  Eq,  (1.  15),  Eq.  (1.17)  may  then  be  optimized 
with  respect  lo  the  power-plant  mass.  This  may  be  done  by  simply  differentiating 
Eq.  (1. 17)  with  respect  to  powerplant  mass  and  setting  the  results  equal  to  zero. 
The  optimum  mass  distribution  for  a  given  value  of  the  payoff,  J,  is  then  given 
by  Eqs.  (1.  lft)  to  (1.21). 

* 

P  / - ’ 

—  “  V<l*l;2)aj  -  aj  . 

ns 


* 


“7  =  1-2  /<i+k  )  aj'  +  aj 
m 
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1.  6  Control  Variables  and  Parameters 

A  control  variable  is  a  variable  which  may  be  manipulated  so  as  to 
control  the  trajectory.  This  may  be  expressed  symbolically  as  Eq.  (1.22). 


dx 

dt 


f (x.u.t) 


In  this  equation,  the  rate  of  change  of  the  state  vector,  x,  is  a  function  of  the 
state  vector,  x,  the  control  vector,  u,  and  the  time,  t..  This  equation  represents 
a  fairly  general  control  problem  and  includes  Eq.  (1. 1)  as  a  special  case.  For  a 
vehicle,  the  control  vector,  u,  is  used  to  control  the  acceleration  of  the  vehicle. 
Examples  of  control  variables  are  throttle  settings,  thrust  magnitudes,  and  thrust 
directions.  The  number  of  control  variables  may  be  greater  or  smaller  than  the 
dimensions  of  the  acceleration  vector.  An  example  of  the  former  is  where  there 
may  be  multiple  engines,  each  with  its  own  control  variables.  An  example  of  the 
latter  is  a  rocket  which  is  controlled  in  three  dimensions  by  a  constant  magnitude 
thrust  at  right  angles  to  the  radius  vector.  Only  one  control  variable,  a  gimbal 
angle,  can  control  all  six  components  of  position  and  velocity  in  an  inverse  square 
field. 

Where  the  number  of  control  variables  is  greater  than  the  dimension  of 
the  acceleration  vector,  it  will  generally  be  possible  to  reduce  the  number  of  con¬ 
trol  variables  by  carrying  out  a  partial  optimization.  This  partial  optimization 
will  maximize  the  acceleration  for  a  given  rate  of  fuel  consumption  or,  equivalently, 


(1.22) 


18 


in 

0 

m 


aj 


The  fuel  mass  and  the  payload  mass  are  then  given  by  Eqs.  (1. 17)  and  (1. 18). 
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(1.16) 


(1.17) 


(1.18) 


These  equations  are  perfectly  general  and  apply  to  rockets  with  either  fixed  or 
variable  thrust,  Equation  (1. 17)  shows  that,  for  a  power-limited  rocket,  the  fuel 
consumption  depends  upon  the  mass  of  the  powerplant.  It  is  possible  for  such 
rockets  to  decrease  fuel  mass  by  increasing  powerplant  mass.  In  general,  there 
will  be  some  optimum  tradeoff  between  powerplant  and  fuel  mass  which  will  maxi¬ 
mize  payload  for  a  given  mission.  For  a  fixed  thrust  rocket,  it  is  necessary  to 
carry  out  the  optimization  of  the  vehicle  design  and  the  trajectory  design  simul¬ 
taneously.  In  this  case,  the  value  of  J  will  depend  upon  the  maximum  acceleration 
available  which  will  depend  upon  the  powerplant  mass.  However,  for  a  variable 
thrust  rocket,  it  is  possible  to  follow  any  required  acceleration  program  with  any 
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minimize  the  rale  of  fuel  consumption  for  a  given  level  of  acceleration.  An  in¬ 
structive  example  of  this  partial  optimization  is  where  a  power-limited  rocket 
and  an  exhaust-velocity-limited  rocket  arc  used  in  parallel.  The  total  mass 
flow  rate  will  be  the  mass  flow  rate  of  the  two  rockets. 


2  2 
•  h  \  *J_ 

'm'  C1  +  2P2  '  C1  +  2P2 


(1.23) 


The  control  variables  will  be  taken  as  the  magnitude  of  the  thrust.  The  thrust  of 
each  engine  will  be  assumed  to  be  variable  from  zero  up  to  some  maximum  value. 


0  <  u  s  f 

1  1  max 


(1.24) 


0  S  u  S  f 

2  2  max 


(1.23) 


f2  max  >  P2/ci 


(1.20) 


The  total  thrust  will  be  the  sum  of  the  thrusts  of  the  individual  engines. 


i  -  v f2  - 


(1.27) 


This  equation  may  be  used  to  eliminate  one  of  the  control  variables  and  to  express 
total  mass  flow  rate  in  terms  of  the  total  thrust  and  the  other  control  variable. 


f-u  u 
2  2 

m  =  ~  +  2P~ 
C1 


(1.28) 


mu  <:nic 


The  mass  flow  rale  as  a  function  of  u^  has  a  single  maximum  whic  h  is  easily  found 
to  be  given  by  Eq.  (1.20), 


* 


(1.20) 


The  exhaust  velocity  at  this  optimum  operating  point  of  the  power-limited-rocket 
is  exactly  twice  the  exhaust  velocity  of  the  exhaust-vclocity-limited  rocket. 


* 


,(1.30) 


With  this  relationship,  it  can  he  seen  that  there  are  three  different  regimes  of 
operation  for  the  two  engines  in  parallel  (sec  Fig.  1.5).  In  the  first  .eghne,  (or 
small  thrusts,  only  the  power-limited  rocket  is  utilized. 


0  e 


(1.31) 


In  the  second  regime,  for  intermediate  thrust,  levels,  the  pov. er-Iiiniicd  rocket  is 
used  at  this  optimum  operating  point  and  the  thrust  of  me  exhaust-velocily-limilod 
rocket  increases  from  zero  to  its  maximum  value. 
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In  the  third  regime,  the  exhaust-velocity-limited  rocket  has  reached  its  maximum 
value  of  thrust  and  the  thrust  of  the  power-limited  rocket  is  increased  to  meet  the 
total  required  thrust. 
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(1.33) 


The  first  and  third  regions  are  examples  of  regions  in  which  the  maximum  occurs 
at  a  boundary  of  one  of  the  control  variables,  in  this  case,  the  thrust  of  the 
exhaust-velocity-limited  rocket.  In  region  two,  we  have  a  stationary  maximum 
for  the  second  control  variable,  the  thrust  of  the  power-limited  rocket. 

There  are  many  possible  formulations  for  the  control  variables  for  a 
given  rocket  propulsion  system.  As  further  examples  of  such  formulations,  some 
of  the  propulsion  systems  of  Section  1.2  will  be  considered.  In  general,  il  will  be 
assumed  that  the  engines  are  gimbaled  so  dial  the  thrust  may  he  pointed  freely  in 
any  direction.  The  control  vector  will  be  taken  to  be  the  thrust  vector  of  the  vehicle. 


I  =  u 


(1.3-1) 


Then  the  mass  flow  rale  may  be  taken  to  be  a  function  of  position,  time,  and  ( lie 
magnitude  of  the  thrust. 


m  -  m  (  r,  t,  u) 


(1.35) 


For  a  chemical  or  nuclear  rocket,  the  mass  flow  rate  will  be  given  by  Kq.  (1.36). 
while  the  control  variable  will  have  a  bound  on  its  absolute  magnitude. 


m 


R 

c 


(1.36) 
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(1.37) 
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If  the  nuclear  rocket  derives  its  source  of  heat  from  isotope  decay  rather  than 
from  a  reactor,  the  decay  of  the  isotope  will  produce  a  continuing  decrease  in 
the  amount  of  power  available,  so  that  the  maximum  value  of  thrust  will  be  given 
by  Eq.  (1.38). 


f 

max 


.  f 0  e-k<l-t0' 


For  an  idealized  nuclear-electric  rocket,  the  mass  flow  rate  will  be 
given  by  Eq.  (1.  39) . 


m  = 


u  •  u 
2P 


In  this  case,  it  will  usually  not  be  necessary  to  set  an  upper  bound  upon  the  magni¬ 
tude  of  the  thrust  as  large  thrusts  will  be  undesirable  because  of  their  exorbitant 
fuel  consumption. 

|  u  |  <■  «> 


If  the  rocket  is  a  solar-electric  rocket  rather  than  a  nuclear-electric  rocket,  the 
available  power  will  be  a  function  of  distance  from  the  sun. 

P  -  P(r) 

For  practical  engines,  where  the  engine  may  only  operate  at  a  single  design  point, 
the  magnitude  of  the  control  vector  may  only  have  two  discreet  values.  In  this 
case,  the  equations  for  the  mass  flow  rate  will  be  the  same  for  a  chemical  rocket, 


(1.38) 


(1.39) 


<1.40) 


(1.41) 
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a  nuclear  rocket  or  n  nuclcar-clectric  rocket. 


=  ft.  f 

max 


A  somewhat  idealized  version  of  a  practical  solar-electric  rocket  may  be  repre¬ 
sented  by  Eqs.  (1.44)  and  (1.45)  . 

|u  | 

m  =  -  J — L 
c 

|u|  s  2211) 

1  C 


(1.42) 

(1.43) 

(1.44) 

(1.45) 


In  this  case,  the  rocket  is  assumed  to  be  capable  of  using  all  the  available  power 
at  constant  exhaust  velocity. 

In  addition  to  control  variables,  the  problem  may  also  contain  a  fixed 
vector  of  control  parameters  which  do  not  change  their  values  wilh  time.  The 
control  problem  is  then  represented  by  Eq.  (1.4G),  where  l  is  the  vector  of 
control  parameters. 

=  T(x,u,/,t)  (1.4G) 


This  problem  arises  when  wc  wish  to  consider  the  optimization  of  the  vehicle  and 
trajectory  for  a  given  space  mission,  such  ns  the  determination  of  the  maximum 
payload  which  can  be  carried  by  a  fixcd-lhnisl  nuclear-elcclric  rocket  on  a  given 
mission.  In  this  case  the  powerplr.nl  size,  which  will  determine  the  value  of  the 
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fixed-thrust  level,  would  be  one  of  the  parameters  of  the  problem.  In  the  same 
example,  the  exhaust  velocity  of  the  rocket  might  be  another  parameter  lobe 
optimized. 


i 


CHAPTER  2:  THE  MAXIMUM  PRINCIPLE 

2.  1  Extrema  of  Functions 

The  problem  of  maximizing  or  minimizing  a  function  is  basic  to  all 
problems  in  space  trajectory  optimization.  The  word  "extremum”  is  used  to 
indicate  a  maximum  or  a  minimum  without  specifying  one  or  the  other.  The 
problem  to  be  considered  in  this  section  is  to  maximize  a  scalar  function  y  of 
a  control  vector  u. 

y  =  f(u)  (2.1) 

The  control  vector  u  is  assumed  to  be  contained  in  a  closed  domain  U.  If  y 
is  a  continuous  function  of  u,  then  the  theorem  of  Wcierstrass  guarantees  that 
the  function  y  always  contains  both  a  maximum  and  a  minimum. 

"Every  function  which  is  continuous  in  a  closed  domain,  U,  / 

of  the  variables  possesses  a  largest  and  a  smallest  value 
in  the  interior  or  on  the  boundary  of  the  domain.  " 

If  the  function  is  differentiable  and  the  maximum  or  minimum  occurs 
at  an  interior  point,  then  the  partial  derivatives  of  y  with  respect  to  each  com¬ 
ponent  of  u  must  vanish.  Such  a  point  where  the  gradient  of  y  vanishes  is 
spoken  of  as  a  stationary  point.  While  an  interior  maximum  or  minimum  of  a 
differentiable  function  must  occur  at  a  stationary  point,  a  stationary  point  need 
not  yield  an  extremum.  For  example,  the  stationary  point  may  correspond  to  a 
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saddle  point  where  both  greater  and  smaller  values  of  the  function  occur  in  the 
near  vicinity  of  Ihc  point.  A  one-dimensional  example  of  a  saddle  point  would  be 
a  point  of  inflection  where  the  first  and  second  derivatives  of  y  with  respect  to 
a  scalar  u  would  vanish.  If  y  possesses  second  partial  derivatives,  then  the 
second  partial  derivative  matrix  would  be  indefinite  for  a  saddle  point,  positive 
definite  for  a  minimum,  and  negative  definite  for  a  maximum.  In  general,  these 
extrema  would  only  be  local  extrema;  that  is,  extrema  within  some  subdomain  of 
U  rather  than  with  respect  to  all  of  IJ. 

The  extrema  of  the  function  y  may  also  occur  on  the  boundary  of  U 
or  at  "corners"  where  y  does  not  possess  partial  derivatives  with  respect  to  u. 
The  theorem  of  Weierstrass  also  applies  to  the  boundary  of  U,  which  is  a  lower 
dimensional  subdomain  of  U.  For  example,  if  the  domain  U  is  a  3-dimensional 
polyhedron,  stationary  extrema  must  be  sought  within  the  interior  of  U,  on  the 
faces  of  the  polyhedron,  and  on  the  edges  of  the  polyhedron.  In  order  to  deter¬ 
mine  the  absolute  extrema  of  y  in  such  a  domain  (if  y  has  continuous  first 
partial  derivatives),  all  of  these  stationary  extrema  must  be  considered  and 
compared  with  the  values  of  y  at  the  vertices  of  the  polyhedron.  If  y  is  a 
linear  function  of  U,  there  will  generally  be  no  stationary  extrema  and  the  ex¬ 
trema  must  occur  at  the  vertices  of  the  n-dimensional  polyhedron.  This  is  the 
basis  of  linear  programming. 
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2.2  Tlie  Optimal  Control  Problem 


An  important  class  of  optimal  control  problems  is  to  maximize  the 
final  value  of  the  payoff  variable  subject  to  the  differential  Eqs,  (2.3)  with 
the  control  vector  u  contained  ir.  some  closed  domain  U. 


=  min. 


=  f(x,u,t) 


The  initial  and  final  times  are  specified  as  are  the  initial  and  final  values  of  the 
state  x.  This  problem  is  a  generalization  of  the  elementary  problem  considered 
in  the  previous  section.  The  cost,  or  payoff,  x*  -epends  on  the  complete  time 
history  of  the  values  of  the  control  vector  u,  rather  than  upon  its  value  at  any 
one  time.  The  payoff  variable  is  what  is  called  a  functional  as  it  depends  upon 
a  continuous  sequence  of  values  of  the  vector  u. 

This  problem  is  much  more  difficult  than  the  elementary  problem  con¬ 
sidered  in  Section  2. 1  and  far  less  is  known  about  its  solution.  There  is  no  general 
existence  theorem  corresponding  to  the  theorem  of  Weierstrass  and  the  sufficiency 
theorems  thai  have  been  developed  are  generally  difficult  to  use  and  are  of  limited 
applicability.  There  is  a  satisfactory  theory  concerning  necessary  conditions  for 
this  problem  and  this  theory  is  the  basis  of  the  analysis  in  this  monograph. 
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2.3  The  Maximum  Principle 


The  maximum  principle  is  a  necessary  condition  for  a  maximizing 


solution  of  the  optimal  control  problem  considered  in  the  previous  section.  Any 


solution  which  does  maximize  the  payoff  must  satisfy  the  maximum  principle. 


-However,  solutions  which  satisfy  the  maximum  principle  may  not  be  maximizing. 
A  solution  which  does  satisfy  the  maximum  principle  will  be  referred  to  as  an 
"extremal . 

The  maximum  principle  is  stated  in  terms  of  a  Hamiltonian  defined  by 

Eq.  (2.4). 


H  =  H (x,  X,  u,  t)  2  X  ■  f(x.  u,  t) 


The  vector  X  is  a  vector  of  Lagrange  multipliers  of  the  same  dimension  as  the 
state  vector  x.  There  is  one  component  of  this  vector  for  each  component  of  the 
state  vector.  The  differential  equations  of  motion  and  the  differential  equations 
governing  the  Lagrange  multipliers  arc  given  by  the  canonical  Kqs.  (2.5)  and  (2.6). 


dx. 


bH 


dt 


<3X. 


i  =  0,  1,  .  . . ,  n 


dX 


bll 


i  =  0,  1,  .  . . ,  n 


dt 


ox. 

l 


These  equations  are  often  stated  in  the  equivalent  vector  form  (2.7)  and  (2.8). 


(2.4) 


(2.5) 


(2.6) 
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dx  _  an 
dt  ~  hi 

dx  SH 

dt  ax 

The  maximum  principle  of  Pontryngin  states  that  the  Hamiltonian  H  of 
Eq.  (2.4)  must  be  maximized  with  respect  to  the  value  of  the  control  vector  u 
contained  within  the  closed  domain  U  at  all  times  from  t°  to  l\  The  maximum 
principle  reduces  the  optimal  control  problem  to  an  infinite  sequence  of  maximiza¬ 
tion  problems  in  Euclidean  space  and  to  a  2-point  boundary  value  problem.  The 
maximization  problem  is  the  problem  of  maximizing  H  with  respect  to  u  at  each 
time.  The  2-point  boundary  value  problem  is  the  problem  cf  determining  the  initial 
values  of  the  Lagrange  multipliers  X.  that  will  cause  the  extremal  to  go  to  the 
correct  final  state. 

Equations  (2.f>)  arc  a  set  of  linear  homogeneous  equations  in  the  Lagrange 
multipliers  X.-  These  equations  are  adjoint  to  the  linearized  variational  equations 
of  Eqs,  (2,5).  A  natural  scaling  of  this  homogeneous  system  of  equations  will  be 
chosen  by  taking  the  terminal  value  of  the  Lagrange  multiplier  associated  with  the 
cost  to  be  unity  (Eq.  (2.9)  ). 


Because  the  Lagrange  multipliers  are  adjoints  to  the  linearized  equations  of  motion, 
and  because  of  the  scaling  used  in  Eq.  (2.9),  they  may  be  interpreted  as  influence 
functions  for  the  cost  and  obey  the  relationship  of  Eq.  (2.10)  at  all  times. 


(2.7) 

(2.8) 


(2.9) 
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The  Hamiltonian  itself  is  the  negative  of  the  influence  function  for  time  and  obeys 


Eg.  (2.11). 


The  system  of  equations  of  Eqs.  (2.5)  and  (2.G)  also  possesses  a  first  integral 
which  gives  the  value  of  the  Hamiltonian  at  all  times  (Eq.  (2. 12)  ). 


H(t) 


H  - 


r  ill 

J  at 

t 


dt 


If  the  system  is  autonomous,  i.e. ,  the  functions  f.  are  independent  of  time,  then 
this  integral  will  be  a  constant  of  the  motion. 

If  all  the  maxima  of  H  with  respect  to  u  are  stationary  maxima,  then 
the  optimal  control  problem  is  reduced  to  a  problem  in  the  classical  calculus  of 
variations.  This  classical  theory  possesses  a  more  satisfactory  and  useful  theory 
of  sufficiency  conditions  than  has  yet  been  developed  for  the  more  complex  optimal 
control  problem. 


2.4  Control  Parameters  and  Free  Boundary  Conditions 

The  problem  considered  in  Section  2.3  may  be  generalized  somewhat  by 
assuming  that  the  boundary  conditions  are  not  fixed  but  may  vary  over  some  ter¬ 


minal  hvpersurface.  This  problem  may  also  be  generalized  by  assuming  that  f 


may  depend  upon  a  constant  vector  of  control  parameters  Z  .  This  makes  the 
Hamiltonian  a  function  of  four  vectors  and  time. 

H  =  X  •  f(x,  u,  Z,  t)  =  H  <X  x.  ",  Z,  t)  (2.13) 


For  any  given  set  of  values  for  the  control  vector  Z,  the  cost  may  be  computed  by 
solving  a  standard  optimal  control  problem.  The  cost  may  then  be  representee'  as 
a  function  of  the  boundary  conditions  and  the  parameter  vector,  Eq.  (2. 14) . 


1  1-0-101- 
x0  =  (x  ,  x  ,  t  ,  t  ,  t) 


(2.14) 


The  problem  of  determining  the  optimal  values  for  the  boundary  conditions  and 
the  parameter  vector  has  now  been  reduced  to  a  problem  in  the  theory  of  extrema 
of  functions.  If  the  cost  given  in  Eq.  (2. 14)  is  differentiable,  then  a  set  of  trans- 
versality  conditions  may  be  derived  which  determine  the  optimal  values  for  the 
boundary  conditions.  For  example,  if  the  initial  conditions  are  fixed  and  the  ter¬ 
minal  conditions  arc  variable,  then  the  variation  in  cost  due  to  variations  in  the 
terminal  conditions  must  be  stationary  and  will  be  given  by  Eq.  (2.  15). 
n 

=  -  Z  x}  ex}  +  II1 6 11  =  0 
i=l 


(2.15) 


This  equation  must  be  satisfied  for  all  variations  which  are  consistent  with  the 
terminal  conditions,  For  example,  if  one  component  of  the  terminal  state  is  un¬ 


specified,  then  by  Fq.  (2. 15)  the  corresponding  component  of  the  Lagrange 


•  •  <  ■ 


multiplier  must  be  zero.  If  the  initial  or  final  time  is  unspecified,  then  the 
Hamiltonian  II  must  be  zero  at  that  time.  More  generally,  if  the  terminal 
state  must  lie  on  some  hypersurface  in  the  space  of  the  state  variables,  then 
the  terminal  Lagrange  multiplier  vector  must  be  normal  to  the  hvpersurface. 
Equation  (2. 15)  may  also  be  used  where  the  terminal  hypersurface  is  a  specified 
function  of  time  and  the  state. 


If  the  payoff  x^  is  differentiable  with  respect  to  the  parameter  vector 

■f,,  then  the  following  conditions  must  hold  with  respect  to  each  component  £  of 

K 

the  parameter  vector. 


r  liL 


dt  =  0 


(2.1G) 


This  condition  may  be  derived  by  considering  each  component  of  the  parameter 
vector  as  a  new  state  variable  whose  derivative  is  equal  to  zero. 


2.5  Singular  Arcs 


The  maximum  principle  determines  the  optimal  control  only  if  the 
Hamiltonian  H  has  a  unique  maximum  value  with  respect  to  the  control  vector  u. 
When  H  has  two  or  more  equal  maxima,  the  problem  becomes  more  complicated. 
One  important  case  arises  when  H  is  linear  in  one  or  more  components  of  u. 

If  the  coefficient  of  the  linear  term  is  nonzero,  the  maximum  of  H  will  occur 
at  the  boundary  of  u.  However,  if  the  coefficient  is  zero,  then  all  values  of  that 


component  of  u  yield  the  same  value  of  H  and  the  maximum  principle  does  not 
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determine  the  control.  This  is  the  ca.se  in  which  singular  arcs  may  arise.  A 


singular  arc  is  a  segment  of  an  optimal  trajectory  of  finite  duration  where  the 
control  does  not  lie  on  the  boundary  and  the  coefficient  of  the  linear  term  in  the 
Hamiltonian  is  identically  zero  over  the  finite  time  interval. 

If  the  problem  is  linear  in  the  state  and  control,  then  the  singular  arc 
will  arise  when  the  solution  is  non-unique  and  may  be  replaced  by  non-singular 
arcs  liaving  the  same  cost.  However,  if  the  problem  is  nonlinear  in  the  state 
x  but  linear  in  one  component  of  the  control  u,  then  singular  arcs  may  repre¬ 
sent  nil  or  part  of  a  unique  minimizing  solution. 

In  practice  it  may  not  be  possible  to  operate  n  control  between  its  maxi¬ 
mum  and  minimum  limits.  However,  an  approximation  to  such  intermediate  con¬ 
trols  may  be  obtained  by  operating  the  control  alternately  at  its  maximum  and 
minimum  values.  If  this  is  done  rapidly,  a  chattering  approximation  to  the  inter¬ 
mediate  control  is  obtained.  Anytime  that  a  Hamiltonian  possesses  two  equal 
maxima,  chattering  may  be  used  to  connect  the  two  maxima  and  the  possibility 
of  a  singular  arc  arises. 

2.  fi  Impulsive  Controls 

If  the  Hamiltonian  is  linear  in  a  component  of  u  and  this  component 
of  u  is  unbounded,  then  a  positive  coefficient  for  this  component  in  the  Hamil¬ 
tonian  implies  that  the  control  should  be  infinite.  The  usual  derivations  of  the 
maximum  principle  do  not  hold  for  this  case.  However,  the  scope  of  the  maximum 
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principle  has  been  extended  so  that  this  important  case  may  he  treated.  The 
details  of  the  results  will  be  presented  in  the  next  chapter  for  the  rochet  problem. 

2.7  Sufficiency 

The  maximum  principle  is  merely  a  necessary  condition  for  an  extremum 
and  its  satisfaction  is  no  guarantee  that  the  solution  to  a  particular  problem  lias  been 
found.  In  special  cases  it  may  be  possible  to  establish  that  ihere  is  a  unique  solu¬ 
tion  that  satisfies  the  maximum  principle.  In  these  cases  the  maximum  principle 
will  be  sufficient  for  an  extremum.  However,  in  general,  it  is  not  possible  to  show 
that  the  problem  possesses  a  unique  solution.  In  fact,  it  is  often  impossible  to 
show  that  the  problem  possesses  any  well-behaved  solution. 

The  only  way  by  which  one  can  determine  absolute  extrema  is  by  calcu¬ 
lating  all  extremals  and  comparing  their  costs.  There  are  a  few  techniques  which 
can  sometimes  be  used  to  determine  local  extrema.  If  the  control  never  lies  on 
a  boundary,  then  the  problem  may  be  treated  by  the  classical  calculus  of  varia¬ 
tions  which  possesses  a  local  sufficiency'  theory  based  on  the  second  variation. 

This  theory  is  a  generalization  of  the  local  sufficiency  theory  for  the  extrema  of 
functions.  If  the  controls  do  lie  on  the  boundary  for  part  of  the  solution,  a  more 
general  approach  nil) si  be  used.  Such  an  approach  ..<js  recently  been  developed 
by  Boltyauski  for  determining  absolute  extrema  within  a  gi.cn  domain.  Unfor  ¬ 
tunately.  this  theory  requires  the  generation  of  th  •  complct:  families  of  extremals 
in  this  domain.  This  is  usually  impractical  if  the  state  vector  x  has  more  than 
two  components  in  addition  to  the  payoff. 
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CHAPTER  3:  GENERAL  THEORY  OF  OPTIMUM 


ROCKET  TRAJECTORIES 

3. 1  Application  of  the  Maximum  Principle 

The  equations  of  motion  of  a  rocket  moving  in  an  arbitrary  time- 
varying  gravitational  field  are  given  by  Eq.  (1. 1)  of  Chapter  1,  This  equation 

will  be  re-written  as  two  first-order  equations,  Eq.  (3.1)  and  (3.2),  so  that  - 

they  will  be  the  first-order  form  to  which  the  maximum  principle  is  applicable. 

v  «■£  +  i  (i",  t)  '  (3.1) 

r  =  V  _  13.2) 

The  rate  of  loss  of  mass  of  the  rocket  is  assumed  to  be  a  general  function  of 
position,  time,  and  the  thrust  vector. 

lia  -  hi  (r,t,f)  (3.3) 

This  system  of  vector  and  scalar  equations  are  the  state  equations  for  a  general 
rocket  traject  ory  optimization  problem.  They  will  form  a  fifth-order  system  for 
planar  flight  in  two  spatial  dimensions  and  a  seven  th-ordcr  system  for  flight  in 
three  spatial  dimensions. 

The  Hamiltonian  of  this  system  is  given  by  Eq.  (3 .  -4) .  This  Hamiltonian 
represents  a  Mayer  formulation  of  the  optimization  problem.  The  vector  \  is 
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nrijoinl  In  the  velocity  vector,  the  vector  ^  is  adjoint  to  the  position  vector,  nnd 
the  scalar  a  is  adjoint  to  the  scalar  mass  of  the  rocket. 


H  =  +  X  •  g  *  Jj  •  V  +  ffiti  (3.4) 

H  -  Il{V,r,m.  X,  ju.  a,  f  ,t)  (3.  5) 

The  Hamiltonian  H  is  a  function  of  the  state  variables,  the  adjoint  variables,  the 
control  vector,  and  the  time.  The  differential  equations  for  the  adjoint  variables 
are  given. by  Kqs.  (3.(1),  (3.7)  and  (3.8).  . 


X  =- 


■Ml 
o  V 


*  U 


<3.(>) 


dr 


o  m 
dr 


(3.7) 


111  -  X  '  f 

dm  2 

m 


<3.  K) 


Because  of  the  simple  relationship  between  the  adjoint  vectors  for  position  and 
velocity.  Kqs.  (3 .  f»>  and  (3.7)  arc  often  written  in  the  second-order  form  given 
by  Eq.  (3.9). 


or 


b  m 

d  r 


(3.9) 


The  adjoint  vecolr  X  plays  an  important  part  in  the  theory  of  optimum  space 
trajectories  and  is  often  called  the  primer  vector,  a  name  originally  introduced 
bv  Lav. den. 
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The  maximum  principle  states  that  the  Hamiltonian,  Eq.  (3.4),  should 
be  maximized  with  respect  to  the  control  vector,  in  this  case  the  thrust  vector. 
Only  two  of  the  terms  in  the  Hamiltonian  are  affected  by  the  thrust  so  that  the 
maximum  principle  reduces  to  the  maximization  of  these  two  terms  (Eq.  (3.10)  ). 


X  *  f  . 

-  +  q  m  -  max 

m 


(3.10) 


Equation  (3. 10)  is  a  general  equation  applying  to  a  great  many  types  of 
rocket  propulsion  devices.  Among  these  types  are  those  where  the  thrust  magni¬ 
tude  must  lie  within  a  given  time-varying  set,  and  to  the  case  where  the  thrust 
magnitude  is  a  fixed  function  of  time,  but  the  direction  may  be  freely  varied. 

/ 

The  one  important  case  not  covered  by  Eq.  (3.10)  is  where  the  maximum  thrust 
magnitude  is  a  function  of  both  position  and  time.  This  case  can  be  conveniently 
handled  by  expressing  the  thrust  as  the  product  of  a  throttle  variable  k,  which 
may  vary  from  zero  to  one,  and  the  maximum  thrust,  of  the  rocket. 


As  the  thrust  is  now  a  function  of  the  position  vector  r,  the  equation  for  the 
primer  vector,  Eq.  (3.9),  must  be  modified  as  Eq.  (3.12). 


X  • 


d K  ,  Am  k 
dr  dr  m 


A 


d  f 

max 

5? 


(3.12) 


At  this  point  it  will  be  assumed  that  the  engine  may  be  rotated  freely 
in  space  so  that  the  maximum  thrust  is  not  a  function  of  direction  and  so  that 
the  fuel  flow  rate  depends  only  upon  the  magnitude  of  the  thrust.  In  this  case, 
Eq.  (3. 10)  implies  that  the  thrust  should  be  aligned  with  the  primer  vector  to 
produce  Eq.  (3.13). 


\f  ^  . 

—  +  am  -  max 

m 


(3.13) 


This  equation  may  now  be  used  to  maximize  the  Hamiltonian  I!  with  respect  to 
thrust  magnitude  as  well  as  thi’ust  direction.  It  is  convenient  to  rc-write  Eq.  (3. 13) 


as  Eq.  (3, 14). 


f 


H  am. 
X 


hi  -  max 


(3.14) 


There  arc  two  important  cases  that  arise  depending  upon  the  curvature  of  the  curve 
of  thrust  versus  fuel  flow  rate.  The  two  cases  are  illustrated  in  Figs.  3. 1  and  3.2. 


In  the  first  case  illustrated  in  Fig.  3.1  the  maximization  of  H  will  generally 
lead  to  a  stationary  maximum  point  away  from  the  boundaries  on  the  thrust. 

This  case  is  characteristic  of  the  idealized  power-limited  rocket.  The  second 
case  illustrated  in  Fig.  3.2  always  has  a  maxima  of  the  Hamiltonian  occurring 
at  the  maximum  and  minimum  values  of  the  thrust.  The  thrust  is  never  operated 
at  intermediate  values  except  in  the  particular  case  where  the  thrust  is  a  linear 
function  of  the  mass  flow  rate.  It  is  possible  to  approximate  to  a  linear  relation¬ 
ship  between  thrust  and  mass  flow  rate  for  these  engines  by  operating  the  engine 
alternately  for  short  periods  of  time  at  both  maximum  and  minimum  thrust,  a 
process  known  as  chattering.  If  the  optimal  trajectory  contains  a  singular  arc, 
such  chattering  behavior  will  form  part  of  the  optimal  solution. 


3.2  Constant  Fxhaust  Vcl.icitv 

In  the  particular  case  of  a  rocket  with  a  constant  exhaust  \clocily,  the 

mass  flow  rale  is  given  by  Hq,  (3. 15),  where  c  is  a  constant,  liquation  (3.  1 5 )  is 

V  " 

then  replaced  by  Fq.  (3. 1C). 


m  - 


f 

c 


f 


g  m 
X  c 


max 


ir  this  case  7  .■  criteria  for  the  magnitude  of  tin-  thrust  are  given  by  Kqs.  (3.  17). 
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X  > 


c  m 
o 


* 

f  =  f 

max 


I 


X  < 


g  m 
c 


* 

f 


=  f 


min 


(o.  17) 


f  .  sfif 
min  max 

If  the  magnitude  of  the  primer  vector  is  greater  than  am/c,  the  thrust  is  turned 
-on  full  throttle.  If  the  magnitude  of  the  primer  vector  is  less  than  jm/c,  the 
engine  will  operate  at  its  minimum  thrust  level,  which  will  lie  taken  to  be  zero. 
“Finally,  in  the  particular  case  of  a  singular  arc  the  magnitude  of  the  primer  vector 
will  remain  equal  to  am/c  along  the  arc,  and  the  value  of  the  thrust  will  lie  be¬ 
tween  its  limits. 


X  * 


a  m 


3.3  Impulsive  Controls 

It  is  interesting  to  determine  what  happens  to  I’qs.  (3.17)  ii  impulses 
are  allowed.  Using  the  results  of  Chapter  1,  Eq.  (1,7),  Eq.  (3.4)  may  be  re¬ 
written  as  Eq.  (3.  18). 

H  •  —  (X  -  —  )  +  X  •  g  -  X  '  V  (3.18) 

me  ' 


The  rale  of  change  of  the  quantity  am  may  be  determined  from  Eq.  (3.19). 
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d(qm) 
d  AV 


mda 


a  dm 


± 

m 


dt 


(3.19) 


When  X  is  less  than  erm/c,  the  thrust  will  be  zero  and  the  rale  of 
change  of  AV  will  be  zero  so  that  there  will  be  no  change  in  the  quantity  a  m  . 
When  X  is  equal  to  erm/c,  the  rate  of  change  of  erm  with  respect  to  AV, 
given  by  Eq.  (3. 19),  will  be  zero.  During  an  impulsive  thrust,  the  value  of  the 
primer  vector  and  its  first  two  time  derivatives  will  not  change.  As  a  result,  the 
right-hand  side  of  Eq.  (3.19)  will  remain  equal  to  zero  during  an  impulsive  thrust 
and  along  a  singular  arc.  This  means  that  as  long  as  the  value  of  thrust  is  allowed 
to  vary  between  zero  and  infinity,  the  quantity  <ym/c  will  be  a  constant  of  the 
motion,  which  may  be  taken  as  unity. 

The  magnitude  of  the  primer  vector  must  be  less  than  unity  when  the 
rocket  is  turned  off  and  must  be  equal  to  unity  when  the  vehicle  is  powered.  It 
can  never  be  greater  than  unity  for  an  optimum  trajectory.  Since  the  primer 
vector  and  its  first  two  time  derivatives  arc  continuous,  all  impulses  between 
the  initial  and  final  times  must  occur  at  local  maxima  of  the  primer  vector.  In 
this  case  the  primer  vector  and  its  derivative  will  be  orthogonal  and  the  last  term 
in  Eq.  (3.  IS)  will  not  change  across  the  impulse.  The  first  term  will  be  zero  on 
both  sides  of  an  interior  impulse  and  may  be  defined  to  be  zero  during  such  an 
impulse.  As  a  result  the  Hamiltonian  for  the  constant  exhaust  velocity  problem 
with  unconstrained  thrust  may  be  taken  as  in  Eq.  (3.20). 


<12 


H  X  *  g  -  A  •  V 


(3.20) 


This  equation  will  not  apply  during  the  terminal  impulses  where  the  Hamiltonian 
will  not  be  defined.  It  will  apply  in  the  open  interval  from  the  initial  to  the  ter¬ 
minal  time. 

3.4  Power-Limited  Rocket 

In  the  case  of  a  power-limited  rocket  where  the  power  may  be  a  function 
of  time  and  position,  the  mass  flow  rate  will  be  given  by  Eq.  (3.21). 

f2 

m  -  -  —  (3.21) 


Equation  {'■).  14)  in  this  case  becomes  Eq.  (3.22). 


o  in 

2  XP 


max 


(3.22) 


This  equation  is  of  the  tjpe  illustrated  in  Fig.  3.  1  and  always  has  the  stationary 
maximum  given  by  Eq.  (3.23). 


* 

f 


X  V 
a  m 


(3.23) 


The  equation  for  the  Lagrange  multiplier  for  the  mass,  a,  is  given  by  Eq.  (3.21) 
and  the  equation  for  the  rate  of  change  of  mass  is  given  by  Eq.  (3.25). 
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A 


Cr  m 


(3.24) 


ill 


P 


2 

m 


(3.25) 


These  equations  may  be  combined  to  yield  Eq.  (3.2G),  which  is  readily  integrated 
to  yield  Eq.  (3.27). 

jT  _  2m 
a  "  r.i 

2  0,  0  2 

erm  =  a  (m  ) 


(3. 20) 


(3.27) 


The  acceleration  of  the  rocket  is  then  given  by  Eq.  (3.26)  and  the  Hamiltonian  by 
Eq.  (3.29). 


* 


c  ill 


(3.26) 


O  m 


X  •  g  -  x  •  v 


2  m 
Hi 


(3.29) 


As  this  is  a  homogeneous  system,  we  may  set  the  Lagrange  multiplier  a  equal 
to  the  value  given  by  Eq.  (3.30)  to  pro: ace  the  Hamiltonian  given  by  Eq.  (3.31). 


,  0  2 

(m  ) 


2  P  - 

H  -  X  —Q  +  X  •  g 
P 


-  X  ■  v  - 


d  ( P  /m) 
dt 


(3.30) 


(3.51) 
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The  quantity  on  the  right  may  be  recognized  as  equal  to  the  quantities  given  by 
Eq.  (3.32).  This  produces  the  final  form  of  the  Hamiltonian  H  for  the  power- 


limited  problem,  Eq,  (3.33). 
d  P°/in 

dt  ~  2  P  m  dt  "  2  p0 

-  -  -  -  x2  P 

11  =  x  • « • x  * v  +  2  7 


It  should  be  noted  that  in  both  cases  where  there  arc  no  bounds  on  the 
thrust,  the  mass  has  disappeared  from  the  problem.  This  is  because  for  these 
problems  the  thrust  may  be  replaced  by  the  thrust  acceleration  as  a  control  vari¬ 
able.  It  should  also  be  noted  that  Eq.  (3.33)  applies  to  an  arbitrary  time-varying 
gravity  field  and  to  a  power  level  that  is  an  arbitrary  function  of  time  and  position. 


3.  5  Summary  of  Results 

From  this  point  on  this  monograph  will  be  primarily  concerned  with  three 
different  trajectory  optimization  problems.  3  he  first  problem,  to  be  referred  to 
as  Problem  Cl,  is  concerned  with  minimum  fuel  trajectories  for  a  constant  ex¬ 
haust  velocity  rochet  having  a  fixed  bound  cn  the  maximum  thrust  level.  In  general, 
the  thrust  level  will  be  assumed  to  be  variable,  but  in  the  special  case  of  minimum 
time  trajectories  the  thrust  will  be  assumed  to  be  fixed  at  its  maximum  value.  In 
such  a  ease  the  minimum  fuel  trajectory  is  also  a  minimum  time  trajectory.  The 
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second  problem,  to  be  referred  to  as  Problem  C2,  will  be  concerned  with  a  con¬ 
stant  exhaust  velocity  rocket  with  no  upper  bound  on  the  thrust  magnitude.  Here 
the  solution  will  consist  of  impulses  and  coasting  arcs  except  in  the  exceptional 
case  of  singular  arcs.  The  payoff  to  be  maximized  will  be  the  negative  of  the 
total  velocity  increment  as  the  velocitj'  increment  should  be  minimized.  For  this 
problem  minimum  time  trajectories  are  meaningless  as  they  require  a  consump¬ 
tion  of  infinite  amounts  of  fuel. 

The  third  problem,  to  be  referred  to  as  Problem  PI,  is  the  minimum 
fuel  problem  for  a  constant  power  rocket  with  unbounded  variable  thrust.  In 
general,  the  power  may  be  a  function  of  position  and  time,  although  most  results 
will  be  for  the  constant  power  case.  For  this  problem  both  minimum  time  and 
time-open  problems  arc  not  of  interest  because  they  correspond  respectively  to 
zero  time,  infinite  fuel  consumption  and  infinite  time,  zero  fuel  consumption 
solutions.  Table  3. 1  summarizes  some  of  the  pertinent  characteristics  of  these 
three  problems. 
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TABLE  3.1 


Problem 

Cl 

C2 

PI 

Exhaust  velocity  -  c 

constant 

constant 

variable 

Maximum  thrust  level  -  f 

max 

bounded 

unbounded 

unbounded 

Power  level  -  P 

variable 

variable 

function  of 
position  and  time 

Payoff  -  x1 

1 

m 

-  AV 

-J 

Terminal  value  of  ^ 

c/m1 

X  '  S  ~  X  '  v 

,  2 

adjoint  to  mass  -  cr 

Hamiltonian  -  H 

1 

x  g-i  ■  v 

/'IN 

ml  .  2 
x-5-x-v+V  - 

m  v  c 

1  p 

Hamiltonian  -  H 


3.6  Boundary  Conditions 


This  section  will  consider  wlwt  effect  various  types  of  boundary  con¬ 
ditions  have  on  the  terminal  values  of  the  Lagrange  multipliers  for  the  three 
problems  defined  in  the  previous  section.  These  boundary  conditions  will  be 
determined  from  the  general  trnnsversalily  condition,  Eq.  (2.  Id)  of  Chapter  2. 
Applying  this  transversality  condition  to  the  rochet  trajectory  problems  of  this 
chapter  produces  Eq.  (3.34).  This  equation  must  be  stationary  for  allowable 
variations  in  the  terminal  conditions. 

-  X  '  6V  - X  •  fir  -  H  5  t 

In  the  case  of  interception  trajectories  the:  terminal  velocity  at  inter¬ 
ception  is  unspecified.  Equation  (3.34)  then  indicates  that  at  interception  all 
three  components  of  the  primer  vector  must  be  zero  (Eq.  (3.  35)  ), 


For  time  fixed  interception  litis  equation  must  be  satisfied  ior  ali  Knee  problems 
Cl,  C2,  and  PI.  Because  the  primer  vector  is  zero  at  the  time  of  interception, 
the  thrust  will  be  zero  at  that  lime  in  all  cases,  l'or  time  open  interception  the 
terminal  position  will  be  a  specified  function  of  time  (Eq.  (3.36)  ). 


The  variation  in  position  duo  to  a  variation  in  interception  time  will  then  be 
given  by  Eq.  (3.37). 
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The  problem  cf  time  fixed  rendezvous  provides  no  freedom  in  the  ter¬ 
minal  conditions  so  that  no  Iransversality  conditions  may  be  imposed.  However, 
there  arc  several  variable  time  rendezvous  problems  of  interest.  It  will  be 
assumed  that  the  target  vehicle  is  falling  free!,,  under  (lie  action  of  gravity. 

Its  position  and  velocity  are  given  by  Eqs.  (3.43),  while  the  variations  in  these 
quantities  arc  given  by  Eqs.  (3.44). 


1  'r4 

!l 

1  U 

v-vt(t) 

(3.43) 

fir  *  Vt(t) 

6  V  -  i  { t ) 

(3.44) 

Applying  Eq.  (3.34)  to  this  problem  yields  Eq.  (3.45). 


d>: 


1 
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(X  *  B 


Vt " »>6t 


(3.45) 


For  the  time  open  case  ior  Problem  Cl  the  stationary  condition  of  Eq.  (3.  15) 
results  in  Eq.  (3.40). 


X 


fr  m 
c 


(3.40) 


This  equation  indicates  that  the  thrust  should  cut  off  at  the  time  oi  rendezvous 
and  that  a  coasting  arc  should  begin  at  lhat  time.  Applying  Eq.  (3.4b)  tr.  Problem 
C2  requires  some  care  because  the  Hamiltonian  may  not  be  defined  at  a  terminal 
impulse.  However,  Eq.  (3.45)  may  be  used  with  the  value  of  the  Ilanul:  >nian 
immediately  prior  to  llie  impulse,  result:  ng  in  Eq.  (3.  17). 
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(3,  '17) 


X  •  (Vt  -  V)  -  X  •  Av  o 

This  equation  may  be  replaced  by  Eq.  (3.48)  as  ahe  primer  vector  must  point 
in  the  direction  of  the  terminal  impulse. 

X  •  X  =  o  (3.48) 

This  equation  indicates  that  the  primer  vector  magnitude  must  have  a  stationary 
maximum  at  the  teiminal  lime.  This  in. plies  tint  for  this  problem  as  well  as 
Problem  Cl  a  coasting  arc  must  begin  at  the  time  of  rendezvous. 

The  minimum  time  rendezvous  problem  is  only  of  interest  for  Problem  Cl. 

Applying  Eq.  (3.45)  yields  Eq.  (3.49', 


This  equation  is  not  very  useful  as  all  it  does  is  determine  a  natural  scaling  few 
the  values  of  the  primer  vector. 

An  interesting  rendezvous  problem  occurs  when  1  he  transfer  time  from 
one  body  in  free  fall  to  another  body  in  free  fall  is  specified,  but  the  time  at  which 
the  maneuver  starts  is  not  specified.  This  case  is  of  interest  for  all  three  problems  — 
Cl,  C2,  and  PI.  Applying  Eq.  (3.45)  to  these  problems  yields,  respectively,  Eqs. 
(3.50),  (3.51),  and  (3.52). 
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(3.50) 


C2 


X  •  AV' 


(3.51) 


PI 


(3.52) 


Many  of  the  interesting  problems  in  space  trajectory  optimization  arc 
problems  of  orbit  transfer.  In  order  to  define  such  problems,  it  is  usually  neces¬ 
sary  to  assume  that  the  gravitational  field  is  only  a  function  of  position  and  not  of 
time.  This  assumption  will  be  made  here.  Such  an  assumption  will  result  in  the 
Hamiltonian  being  constant  for  Problems  Cl,  C2,  and  those  examples  of  Problem  PI. 
where  the  power  is  independent  of  the  time.  The  variations  in  position  and  velocity 
on  the  target  orbit  will  be  given  by  Eqs.  (3.53),  where  the  variable  r  refers  to  a 
fictitious  motion  on  the  target  orbit  and  does  not  correspond  to  time  on  the  actual 
trajectory. 

6  r  =  Vt6r  6  V  =  g6r  '  (3.53) 

For  time  fixed  orbit  transfer  for  Problem  Cl,  Eq.  (3.34)  yields  Eq.  (3.54)  which 
will  be  true  at  both  ends  of  the  trajectory. 


X  •  g  -  X  *  V  =  0 


(3.54) 
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Because  die  Hamiltonian  will  be  constant  for  this  problem,  Eq.  (3.55)  may  also 
be  obtained. 


"  -  £  b 


ni 


t 
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For  Problem  C2  with  time  fixed  the  variations  in  position  and  velocity  will  be 
assumed  to  be  given  by  l'tjs.  (3. 50) . 


6  r  «  (V  i  AV)6r 


6  V  --  g  6  r 


This  results  in  Kqs,  (3.57)  and  (3.58). 


fix.  «  0  =  (X  *  g  -  A  •  V  T  X  '  AV)  fir 


11  -  -  A  *  AV 


=  +  A  *  AV 


F or  Problem  PI,  Eq.  (3.  hi)  is  again  obtained  and  if  the  power  is  not  a  function 
of  time,  Kq.  (3.59)  may  also  be  obtained. 


H  =  -r- 


0 


l  JL  ■*' 
0 


2 


For  time  open  orbit  transfer  for  Problem  Cl,  Eq.  (3.54)  still  holds 


and,  in  addition,  Eq.  (3.  GO)  holds.  The  Hamiltonian  will  be  zero. 


(3.  GO) 


The  terminal  times  for  this  time  open  problem  must  represent  the  beginning  of 
coasting  arcs. 

For  time  open  orbit  transfer  for  Problem  C2,  Eq.  (3.58)  must  hold. 
However,  in  this  case,  because  the  time  is  not  specified,  the  Hamiltonian  must 
be  zero.  Therefore,  Eq.  (3. 61)  must  hold. 


X  •  X  =  0 


(3.61) 


Once  again,  as  in  Problem  Cl,  the  terminal  powered  arcs'  must  join  onto  coasting 


arcs. 


For  minimum  time  orbit  transfer  in  Problem  Cl,  Eq.  (3.54)  must 
hold  and  Eq.  (3. 62)  must  hold.  ^ 


1  in 
A  =  c 


(3.62) 


As  for  minimum  time  rendezvous,  this  provides  a  natural  scaling  for  the,; I.agrange  . . . 

multipliers. 

Another  interesting  problem  is  the  problem  of  reaching  a  given  energy. 

The  energy  in  a  time  varying  gravitational  field  may  be  represented  as  in  Eq.  (3.63) 
where  U  is  the  lime  varying  gravitational  potential. 


V  •  V 

E  =  2~~  +  U  (r,  t) 


(3.63) 


For  the  time  fixed  problem  the  allowable  variations  in  velocity  and  position  must 


satisfy  Eq.  (3.  G-i). 

v  ■  Fv  -  g  •  fi  r  «  o 

The  change  in  cost  from  Eq.  (3.34)  is  given  by  Eq.  (3.03), 

6  xj  -  \  •  rv  -  i  • « ?  =  o 

In  order  for  both  (3.GJ)  and  (3.  Go)  to  be  satisfied,  Eqs.  (3.GG)  must  bo  satisfied. 


For  all  three  problems  the  quantity  on  the  left-hand  side  of  Eq.  (3.  07)  is  equal 
to  the  right-hand  side. 

x  '  i  -  x  •  v  =  i  •  (vt- v) 

For  both  Problem  Cl  and  PI  the  velocity  is  continuous  so  that  both  quantities 
will  be  equal  to  aero.  For  Problem  C2  the  value  of  II  will  not  be  zero  but 
will  be  given  by  the  right-hand  side  of  Eq.  (3.G7).  It  should  be  noted  that  the 
terminal  impulse  for  Problem  C2  will  be  in  the  direction  of  the  velocity  vector 
so  that  V  and  V  will  be  colinenr. 

To  consider  the  problem  of  reaching  a  given  energy  with  variable 
terminal  time,  it  will  be  assumed  that  the  gravitational  potential  is  not  a  function 
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(3.04) 


(3.  G5) 


(3.GG) 


(3.  G7) 


A 


of  lime.  For  the  time  open  problem  of  reaching  a  given  energy  for  Problem  Cl, 
the  trans versalily  conditions  will  result  in  Eqs.  (3.66),  (3. 67)  and  (3. 60).  For 
Problem  C2  the  corresponding  equations  Will  be  (3.66),  (3.67)  and  (3.61).  This 
implies  that  for  this  problem,  Eq.  (3. 68)  must  be  satisfied. 


This  equation  has  a  simple,  physical  interpretation.  The  rate  of  change  of  energy 
is  proportional  to  velocity  so  that  the  terminal  orbit  must  be  entered  at  a  point 
where  the  velocity  is  a  maximum.  The  stationary  condition  for  this  to  be  true 
is  given  by  Eq.  (3.  68),  which  will  be  satisfied  at  an  apse  of  the  terminal  orbit. 
For  minimum  time  transfer  to  a  given  energy  level  in  Problem  Cl,  the  trans- 
versalily  conditions  will  result  in  Eqs.  (3.66),  (3.67)  and  (3.62). 


(3.68) 


CHAPTER  4:  TRAJECTORY  OPTIMIZATION 
IN  FIELD-FREE  SPACE 


4. 1  Integration  of  the  Adi  pint  Equations 


If  the  mass  flow  rate  of  the  rocket  is  independent  of  position,  then  the 


adjoint  equations  become  uncoupled  from  the  equations  of  motion  for  linearized 
gravity  fields.  An  important  special  case  of  such  fields  is  field-free  space  far 
from  any  massive  body. 


In  field-free  space  the  differential  equation  for  the  primer  vector  is 


given  by  Eq.  (4.1) 


X  =  0 


This  equation  is  readily  integrated  to  yield  Eq.  (4.2). 


X  =  X°  +  X°  (l-t°) 


The  tip  ol  the  primer  vector  moves  along  a  straight  line  at  constant  velocity. 
This  line  may  be  spoken  of  as  the  locus  of  the  primer  vector.  The  primer 
vector  always  lies  in  a  plane  determined  by  the  origin  and  tins  locus.  The 
magnitude  of  the  primer  vector  is  given  by  Eq.  (4.3). 


-0  -0-0-0  0  -0-0  02 

X  -x  +  2 x  -X  (t-t  )  +  x  -x  (t-t  ) 
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'.1  his  equation  represents  a  hyperbola  which  is  concave  upwards.  Because  of  (he 
shape  of  this  curve  there  will  be,  at  most,  one  coasting  arc  for  Problem  C2  and 
cither  one  or  two  impulses.  For  Problem  C2  there  will  be,  at  most,  one  coast¬ 
ing  arc  with  either  one  or  two  thrusting  arcs.  If  the  final  position  is  open,  the 
derivative  of  the  primer  vector  will  be  identically  zero  and  the  primer  vector 
will  be  constant  in  magnitude.  This  will  be  a  singular  case  for  Problems  Cl  and 
C2,  the  thrust  direction  will  be  specified  but  the  time  of  application  and  magni¬ 
tude  of  the  thrust  will  be  unspecified. 

4.2  Constant  Exhaust  Velocity  Propulsion  with  Unbounded  Thrust  Magnitude 

i 

For  interception  the  previous  .chapter  has  shown  that  the  primer  vector 
must  lie  zero  at  the  terminal  time,  Eq.  (4.4). 

!  X1!  -  0  (4.4) 

The  primer  vector  can  only  beeonic  zero  if  its  rate  of  change  is  directed 
antiparallel  to  its  initial  direction.  This  result  implies  that  the  direction  of 
thrust  during  interception  must  be  constant.  This  constant  thrust  direction  for 
interception  is  independent  of  any  consideration  of  the  type  of  propulsion  and  holds 
for  all  propulsion  systems  in  field-free  space. 

The  magnitude  of  the  primer  vector  must  decrease  linearly  to  zero 
at  interception.  There  will  lie  only  one  impulse  represented  by  Eq.  (4.5),  which 
may  be  rc-wriltcn  as  (4.G). 
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-  1  -  0 

AV  “  irr  ■  v° 


-o  -o 

-  rt  "  1  -0  -0 

AV  =  +  vt  -  V 

t  ■*  t 


The  magnitude  of  the  required  velocity  increment  may  be  found  by 
taking  the  dot  product  of  Eq.  (4.  (>)  with  itself.  The  velocity  increment  is  a  func¬ 
tion  of  the  initial  relative  position,  the  relative  velocity,  and  the  total  time 
allowed.  For  given  initial  conditions  there  is  an  optimum  time  for  interception 
which  minimizes  the  velocity  increment.  This  optimum  time  may  be  found  by 
determining  the  stationary  point  of  the  norm  of  Eq.  (4.G).  Carrying  out  this 
computation  yields  the  two  minimizing  roots  given  by  Eq.  (4.  7). 


t!-i° 


o  -0  -0  -0 

( " r  )  *  ( rt  - r  ) 

To  To  To  rr  ’ 

(rt  -  r  )  •  (Yt  -V  ) 


The  left-hand  root  represents  the  optimum  interception  time  when  the  vehicle  is 
approaching  the  target,  while  the  right-hand  root  represents  the  optimum  inter¬ 
ception  time  when  the  vehicle  is  receding  from  the  target.  The  left-hand  root 
lias  an  interesting  physical  interpretation  ns  it  is  the  negative  of  the  range  to  the 
target  divided  by  the  rate  of  change  of  range.  The  minimum  velocity  increment 
for  these  optimum  times  is  given  by  Eq.  (4.8),  the  left-hand  expression  again 
applying  to  vehicles  approaching  the  target  and  the  right-hand  expression  for 


vehicles  receding  from  the  target. 
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AV 


(4.8) 


The  left-hand  term  in  Eq.  (4.8;  is  equal  to  the  magnitude  of  the  relative 
velocity  multiplied  by  the  sine  of  the  angle  between  the  relative  velocity  vector 
and  the  line-of-sight  to  (lie  target.  It  applies  only  when  the  vehicle  is  approaching 
the  target.  When  the  vehicle  is  receding  from  the  target  the  optimum  maneuver 
is  to  stop  the  relative  motion  except,  for  an  infinitesimal  residual  velocity  which 
will  accomplish  interception  as  the  time  becomes  infinite. 

When  the  vehicle  is  moving  towards  the  target,  the  optimum  intercep¬ 
tion  has  several  interesting  geometric  properties.  As  required  by  Eq.  (8.89)  of 
the  previous  chapter,  the  rate  of  change  of  the  primer  vector  must  be  perpendic¬ 
ular  to  the  relative  velocity  at  interception.  The  thrust  is  also  directed  from  the 
unperturbed  position  that  the  interceptor  would  have  at  the  terminal  time  to  the 
position  that  the  target  has  at  the  terminal  time.  This  last  result  will  be  tru6 
for  any  interception  in  field-free  space  with  any  piopulsion  system. 

Another  interesting  interception  problem  is  when  the  transfer  time  is 
fixed  but  the  time  at  which  the  initial  orbit  is  left  is  not  specified.  The  required 
velocity  increment  for  this  ease  is  given  by  Eq.  (4,9),  which  may  be  re-written 
as  Eq.  (4. 1C). 


GO 
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For  rendezvous  there  will  be  two  impulses  which  will  always  be  given 
at  the  beginning  and  end  of  the  maneuver.  The  first  impulse  puts  the  vehicle  on 
a  collision  course  with  the  target.  This  is  the  same  impulse  as  is  required  for 
interception.  The  second  impulse  cancels  the  relative  velocity  between  the  two 
vehicles  at  interception.  The  total  velocity  increment  required  is  given  by 
Eq.  (-l.Fl). 


AV 


- 0  -0  -0  -0  1  0, 
rt  -r  +(Vt  -V  )(l  -t  ) 
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rt  -r  1 

t1 .1° 


(4.13) 


The  magnitude  of  the  second  impulse  depends  only  upon  the  initial  range  to  the 
target  and  the  allowable  rendezvous  time.  The  total  velocity  increment  required 
for  rendezvous  is  a  monoionically  decreasing  function  of  the  time  allowed.  The 
optimum  time  for  rendezvous  is  always  infinite.  This  optimum  rendezvous  time 
is  the  same  as  the  optimum  time  for  interception  when  the  vehicle  is  receding 
from  the  target. 

If  the  transfer  time  is  fixed  but  the  time  of  initiating  the  maneuver  is 
open,  the  rendezvous  may  be  analyzed  by  the  same  type  of  analysis  as  was  used 
for  the  corresponding  interception  case.  Once  again,  the  first  impulse  will  be 
the  same  as  that  required  for  interception  and  the  second  impulse  will  cancel 
the  relative  velocity  difference.  The  total  velocity  increment  is  given  by 


Eq.  (4.14). 
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(4.14) 


The  optimum  time  for  initialing-  the  rendezvous  maneuver  may  be  found  by  dclcr- 
mi.nhig  l.hr  stationary  minimum  of  Eq.  (4.14)  with  respect  to  the  terminal  time. 
The  optimum  terminal  time  is  given  by  Eq.  (4.15). 


-0  -0  -0  -0 
(rt  -  r  )  *  <Vt  -  V  ) 

To  -  —  ~ 

(vt  -v  )  *  (vt  -V  ) 


(4.15) 


Unlike  the  corresponding  interception  case,  the  times  to  start  and  end  the  maneuver 
are  centered  around  the  unperturbed  closest  approach  time  rather  than  terminating 
at  that  time.  The  total  velocity  increment  for  this  case  is  given  by  Eq.  (4. 16). 


AV  V  (Vt°  -  V°)  •  (V®  -  V°)  +  4 
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(4.16) 


The  two  velocity  increments  arc  equal  in  magnitude.  They  each  have  a  component 
normal  to  the  relative  velocity  whose  magnitude  is  the  same  as  that  for  interception 
in  the  same  transfer  lime.  These  two  velocity  components  cancel  each  other. 

Each  velocity  increment  also  has  a  component  which  is  equal  to  half  of  the  relative 
velocity  and  opposed  to  the  relative  velocity  direction. 

This  minimum  velocity  increment  for  rendezvous  with  a  fixed  transfer 
time  is  the  same  as  the  minimum  velocity  increment  for  transfer  from  the  first 
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At  the  terminal  time  the  position  of  the  vehicle  must  be  the  same  as  the  position 
of  the  target.  Using  this  fact  and  substituting  Eq.  (4.20)  into  Eq.  (4.1R)  yields 
Eq.  (4.21). 
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(4.21) 


Equations  (4.21)  and  (4.20)  provide  the  solution  of  the  boundary  value  problem 
and  allow  the  determination  of  all  quantities  of  interest,  including  the  payoff. 
Eq.  (4.22). 
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(4.22) 


This  expression  for  J  has  an  .absolute  minimum  of  zero  at  infinite  time.  It 
also  has  a  local  minimum  with  respect  to  t*  if  the  angle  between  the  realtive 
velocity  and  position  vectors  is  less  than  30°  and  the  intercepting  vehicle  is 
approaching  the  target.  The  time  at  which  this  local  minimum  occurs  is  given 
byKq.  (4.23). 
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During  an  optima]  interception  the  thrust  acceleration  will  decrease 
linearly  to  zero.  However,  for  guidance  purposes  it  is  useful  to  have  a  closed 

I  j 

loop  synthesis  of  the  optimal  control  in  terms  of  the  current  stale  of  the  vehicle, 
the  current  slate  of  the  target,  and  of  the  time  to  go.  Such  a  synthesis  is  easily 
achieved  by  regarding  the  initial  state  in  Eqs.  (4,20)  and  (4.21)  as  the  current 
state  to  yield  Eq.  (4.24). 


(4.24) 


The  numerator  in  Eq.  (1.22)  for  the  payoff  is  simply  the  square  of  the 
unperturbed  distance  between  the  two  vehicles  at  the  terminal  time.  For  inter¬ 
ception  with  a  given  transfer  time  but  a  variable  terminal  time,  the  optimum 
terminal  time  will  obviously  be  the  time  of  closest  approach.  This  is  the  same 
as  the  optimum  terminal  time  for  the  impulsive  ease,  Kq.  (4.11).  The  minimum 
payoff,  J.  for  this  case  will  then  be  given  by  Eq.  (4.25). 


(Vt°-V°)  X  (?t°-  r°)| 
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(4.25) 


For  rendezvous,  Eqs.  (4.17)  and  (4.  IS)  must  be  solved  simultaneously 
to  yield  the  initial  values  of  the  adjoint  vector  given  by  Eqs.  (4. 2(5)  and  (4.27). 
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The  corresponding  value  of  the  payoff  J  is  given  by  Eq.  (4.28). 
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This  payoff  J  is  :i  monotonically  decreasing  function  of  the  transfer  time.  The 
optimum  time  for  rendezvous  is  once  again  infinite  as  it  was  for  the  impulsive 
case. 

The  optimum  control  for  rendezvous  may  be  synthesized  exactly  as  it 
was  for  the  interception  case  to  yield  Eq.  (4.29). 
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In  this  case  the  magnitude  of  the  acceleration  during  an  optimum  rendezvous  will 
vary  hyperbolically  with  time. 

For  the  case  where  the  transfer  time  is  specified  but  the  lime  of  starling 
or  terminating  the  maneuver  is  not,  the  payoff  is  given  by  Eq.  (4.30). 
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This  equation  has  a  stationary  minimum  at  precisely  the  same  time  as  that  for 
impulsive  rendezvous  given  by  Eq.  (4.15).  The  corresponding  payoff  for  this 
case  is  given  by  Eq.  (4.31). 
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(4.31) 


The  first  term  in  Eq.  (4.31)  represents  the  cost  of  a  constant  deceleration  which 
cancels  the  relative  velocity  between  the  two  vehicles.  The  second  term  repre¬ 
sents  the  cost  of  a  change  in  position  with  no  change  in  velocity  at  right  angles 
to  the  velocity  change.  The  latter  thrust  component  varies  linearly  with  time 
and  passes  through  zero  at  the  unperturbed  closest  approach  time. 
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Substituting  Eq.  <5. 11)  into  Eq.  (5.10,.  produces  Eq.  (5.17). 

^a-?+2r-i)  +  3M.(xi..  o^)  =  0 


Eq.  (5. 17,  ,  mny  be  integrated  with  respect  to  lime  to  produce  the  ndditional 
integral  0f  motion  for  these  two  problems _ Eq.  (5.  IS). 


(5.17) 


*'?  +  2 ;-X  +  3,lt  -  (  AV,  5J)  =  const. 


TIki  integration  can  be  carried  out  for  the 


constant  exhaust  velocity  case  be¬ 


cause  the  thrust  is  turned  on  only  when  the  primer  vector  has  unit  magnitude. 

The  integral  holds  on  the  open  interval  from  the  initial  time  to  the  final  time 
even  if  there  are  singular  arcs. 


(5.18) 
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5.2  Variation  of  Parameters 


For  many  problems  it  is  convenient  to  rewrite  the  equations  of  motion 
in  terms  of  Die  variation  of  orbital  elements,  instead  of  the  position  and  velocity 
vectors  used  in  the  previous  section.  The  oi’bital  element  formulation  is  useful 
for  both  coasting  arcs  where  the  vehicle  is  unpowered  and  for  low  thrust  tra¬ 
jectories  where  the  elements  of  the  orbit  may  be  expected  to  vary  slowly.  Eq. 


(f>.2)  is  replaced  by  six  first  order  equations  for ; the  rates  of  change  of  the 
orbital  elements.  The  elements  chosen  ai*e  semimajor  axis  a,  eccentricity  e, 
argument  of  perigee  to,  inclination  i,  longitude  of  the  node  0,  and  mean  anomaly 
M.  Results  for  other  sets  of  elements  can  be  derived  in  the  same  way.  The 
rates  of  change  of  these  elements  may  be  found  in  any  standard  text  celes¬ 
tial  mechanics  and  are  given  by  (5. 19)  through  (5.24). 
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a~l  F  (P  cos  to  -  f  f  sin  05) 
z  4  15  _ 


mf. 
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dM 

dt 


Q&4  IQ'  F  (f  f  uu  -I-  X  £3 ill  w./ 

_  „  /_  1  5 _ ± _ 

dt  vy  mf,  sin  i 
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F  Fr(t4^  '  2e4»  -  *9  W.  +  £ 


+  - 


(5.  231 


(5.24' 


mef 


The  me nn  and  eccentric  anomalies  are  related  by  Kepler's  Eq.  (5.25)  and  the  f. 
arc  given  by  Kqs.  (5.2G). 


M  =  E  -  e  sin  E 


(5.25) 


f1  s  sin  E 


fg  E  cos  E 


f  s  1  -  e  cos  E  ^ 
3 


r  =  cos  E  -  e  > 
4 


(5.2G) 


The  Hamiltonian  for  this  formulation  is  given  by  Eq.  (5.27). 

•  •  • 

H  =  X  a  +  X  e  +  X  to  +  X.i  +  XfP  +  X.  M  +  erm 
a  e  co  i  m  m 

/ 

This  Hamiltonian  is  the  same  as  the  Hamiltonian  given  by  Eq.  (5. 1).  The  trans¬ 
formation  between  the  two  formulations  may  be  found  by  the  standard  techniques 
of  canonical  transformation.  A  different  technique  will  be  utilized  herein.  The 


(5.271 
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basis  of  Dus  technique  is  noting  that  the  instantaneous  values  of  X  and  X  are 
independent  of  the  thrust  level  so  that  X  and  X  have  the  same  value  on  the 
osculating  unpowered  conic  and  on  the  powered  trajectory.  For  the  os¬ 
culating  conic  the  thrust  will  be  zero  which  yields  Eq.  (5.28)  when  Eq.  (5. 1) 
and  (5.27)  arc  equated. 


X-r  3  -  X-r  -  3 


(5. 281 


The  components  of  the  primer  vector  in  the  cylindrical  coordinate 
systems  whose  axis  is  perpendicular  to  the  orbital  plane  may  be  determined  by 
equating  the  components  of  the  primer  vector  in  Eq.  (5. 1)  and  (5.28). 


P'A*  W*2 

"  vy  r 


(5.29) 


(5.30) 


(5.311 


The  quantities  ft  and  ft  are  given  by  Eqs.  (5,32). 

X  C* 


ft,  -  X.  cos  to  4  X«  sin  to  esc  i  -  X  sin  CO  cot  i 
1  1  ft  to 


ftn  -  -X.  sin  to  4  \  cos  to  esc  I  -  X  cos  to  cot  i 
2  i  ft  to 


0.  oi 


J 


The  derivative  of  the  primer  vector  is  calculated  by  determining  its 
value  on  the  osculating  unpowered  ellipse.  The  Hamiltonian  on  this  ellipse  is 
given  by  Eq.  (5.33). 


With  tliis  Hamiltonian,  the  rates  of  change  of  the  Lagrange  multipliers  for  all 
the  elements  except  the  semimajor  axis  ai'c  zero.  The  rate  of  change  of  the 
Lagrange  multiplier  for  the  semimajor  axis  is  given  by  Eq.  (5.34). 


(5.33) 


x  =-\ 

»  2  M/a5 

The  components  of  the  derivative  of  the  primer  vector  in  the  r,  0,  z  directions 


(5.34) 


arc  given  by  Eqs.  (5.35),  (5.3(5),  and  (5.37). 

•  BX  ,  c)X  >  • 

\  =  i#E+srVxe9 
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(5.  35). 


ax 

X0=  ~ 

ax  .  ax  . 

7  7 

X  =  TF  E  +  TT 
z  3  E  oX  a 
a 

The  first  two  terms  on  the  right-hand  side  represent  the  rates  of  change  due  to 
the  variables  in  the  equations  while  the  last  terms  on  the  right-hand  side  are  due 
to  the  rotation  of  the  coordinate  system.  The  rates  of  change  of  the  eccentric 
anomaly  and  central  angle  are  given  by  standard  two-body  formulas. 
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(5.36) 

(5.37) 
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Carrying  out  these  operations  results  in  Eqs.  (5.38),  (5.39),  and  (5.40). 


-2X  ae  -X  ef0f  -  (X  -  f  X  Iff  +  X  e2f  f 
X  =  —  a  g  2  5  i  _5_M  15  AM  \l3 
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-^erifc  +  (X  -  fA„)  fn  -i-  X  ef 
-  e  15  6l!  5  M'  2  M  5 


(5.  38) 


(5.39) 
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“)S  f  f  +3  f 

A.  =  11  0  6  & 

z  •«,«. 

These  six  expressions  for  the  primer  vector  and  its  rates  of  change  are  general 
expressions  which  are  true  both  for  powered  and  unpowered  orbits.  They  repre¬ 
sent  the  results  of  a  canonical  transformation  between  the  two  forms  of  the  equa¬ 
tions  of  motion.  In  the  particular  case  of  an  unpowered  coasting  arc,  the  oscula¬ 
ting  orbit  will  be  identical  with  the  coasting  arc  at  all  points  and  Eq.  (5.34)  may 
be  integrated  to  yield  Eq.  (5.41). 


(5.30) 


vxM 

On  a  coasting  arc,  the  Lagrange  multipliers  for  all  elements  except  the  semi¬ 
major  axis  will  be  constant.  It  should  also  be  noted  that  the  equations  for  the 
primer  vector  and  its  derivative  on  a  coasting  arc  are  the  same  as  the  vari¬ 
ational  equations  for  position  and  velocity  respectively.  This  means  that 
any  of  the  innumerable  solutions  of  the 
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variittion.nl  equations  of  two-body  orbits  may  be  used  to  determine  the  primer 
vector  time  history  during  coast.  The  present  formulation  has  the  possible 
advantage  of  identifying  the  integration  constants  in  terms  of  the  Lagrange 
multipliers  of  the  orbital  elements. 

The  integrals  of  motion  derived  in  section  (5. 1)  will  now  be  expressed 
in  terms  of  the  orbital  elements  and  their  multipliers.  This  will  be  done  by 
direct  calculation  of  the  dot  and  cross  products  of  the  primer  vector  with  the 
position  and  velocity  vectors.  The  components  of  position  and  velocity  are 
given  by  the  following  standard  two-body  equations. 


The  first  constant  of  the  motion,  the  Hamiltonian,  is  given  by  Eq.  (5.42)  for 
the  constant  exhaust  velocity  problem  and  for  the  constant  power  problem. 


It  should  be  noted  that  for  the  constant  exhaust  velocity  problem,  this  provides  a 
determination  of  in  terms  of  the  current  state.  The  value  of  (lie  primer 
vector  for  the  power  limited  problem  may  be  determined  by  squaring  and  adding 
Eqs.  (5.29),  (5.30)  and  (5.31). 

The  vector  constant  will  be  determined  by  direct  calculation  to  be 
given  by  Eq.  (5. 43)  where  the  components  of  the  column  vector  on  the  right 
side  are  in  the  r,  9,  and  /.directions. 
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rxX  +  Xxr  = 


f  R  +  f  f  8 
4P1  15  2 


f  f  8  -  f  8 

1  5P1  4P2 


-(5.43) 


This  result  is  precisely  the  result  obtained  by  Moyer  (A1AA  J.  7,  1232-1235) 
for  the  time-open  constant  exhaust  velocity  problem.  It  is  shown  in  that  paper 
that  Eq.  (5.43)  yields  the  following  equations  for  the  Lagrange  multipliers  of 
the  Euler  angles  of  the  orbit. 


(5.441 


cos  i  +  sin  i  (C^  sin  Cl  -  cos  0) 


X.  =  C„  sin  Cl  +  C,  cos  O 

l  2  1 


The  two  constants  in  these  equations  are  given  by  Eqs.  (5.47)  and  (5.48). 


(5.45> 


(5.4G) 


C1  "  \  cos  Q°  “  ^  cot  *°  sin  Cl°  -i-  X°  sin  Cl0  esc  i° 


(5.471 


C2  “  X.  sin  Cf  +  Xq  cot  i°  cosd°  -X°  cos  fi°  Csc  i 


,0  0 


(5.48) 


It  should  be  noted  that  equations  (5.44),  (5.45),  and  (5.46)  are  simply  the  dot 


products  of  the  vector  constant  with  the  axes  of  rotation  of  the  Euler  angles. 


Finally,  the  quantity  in  the  last  integral  may  be  determined  by  direct 


calculation  to  be  given  by  Eq.  (5.49). 

<  • 

X-  r  +  2r’  X  =  -2\  a 
a 


(5. 49' 


Equation  (5.18)  then  provides  a  direct  representation  of  the  Lagrange  multi¬ 
plier  for  the  semimajor  axis  in  terms  of  the  Hamiltonian,  the  time,  the  cost, 

/ 

and  the  semimajor  axis.  It  should  be  noted  that  Eq.  (5. 18)  reduces  to  Eq.  (5.41) 
on  a  coasting  arc. 
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5.3  Maxima  of  the  Primer  Vector 

/ 

For  the  constant  exhaust  velocity  case,  the  thrust  may  only  be  applied 
at  absolute  maxima  of  the  primer  vector  having  unit  magnitude. 

X-X-l  (5.50) 

Since  the  primer  vector  is  continuous  and  has  continuous  first  and  second  time 
derivatives,  all  maxima  except  those  at  the  initial  and  final  time  must  be  sta¬ 
tionary  maxima. 

X'X=0  (5.51) 

A  necessary  condition  for  these  stationary  values  to  be  local  maxima  is  that  the 
second  derivative  of  Eq.  (5.50)  with  respect  to  time  must  be  negative. 

•  •  ••  ; 

X'X+X-X^O  (5.52) 


The  dot  product  of  the  primer  vector  with  its  second  derivative  may  be  found 
from  Eq.  (5.4)  to  be  given  by  Eq.  (5.53) 

X.'i  =  -£(1.3\a)  (5.53) 

3  r 

r  / 

As  the  dot  product  of  any  vector  with  itself  is  necessarily  positive,  Eq.  (5.52) 
yields  the  additional  inequality  of  Eq.  (5.54). 


i:i<  o 


(5.54) 
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Equations  (5.53)  and  (5.54)  together  yield  Eq.  (5.55). 

X2  =•  sin2  35. 2G4°  (5.55) 

r  o 

A  necessary  condition  for  Die  primer  vector  to  have  a  stationary  maximum  is 
that  the  thrust  direction  must  be  within  35.  2G4  degrees  of  the  local  horizontal 
direction. 

Equation  (5.52)  may  be  used  to  obtain  a  bound  on  the  rate  of  change 
of  the  primer  vector  given  by  Eq.  (5.5G). 

X-X  ^(i-3\r)  (5.5G) 

r 

The  maximum  rate  of  rotation  of  the  primer  vector  at  a  stationary  maxima  is 
equal  to  the  rate  of  rotation  of  a  satellite  in  a  circular  orbit  at  the  same  radial 
distance.  This  maximum  rotation  rate  can  only  be  realized  if  the  radial  com¬ 
ponent  of  the  primer  vector  is  zero.  If  the  radial  component  of  the  primer 
vector  has  the  maximum  allowable  value  given  by  Eq.  (5.55),  then  the  rate 
of  rotation  of  the  primer  vector  must  be  zero.  > 

For  the  case  of  singular  arcs,  the  magnitude  of  the  primer  vector 
must  remain  unity  over  a  finite  time  interval.  In  these  cases,  the  inequality 
in  Eqs.  (5.52)  and  (5.5G)  must  be  replaced  by  an  equality. 
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Chapter  (>:  LINEARIZED  POWER-LIMITED  TRANSFER 
IN  THE  VICINITY  OF  AN  ELLIPTIC  ORBIT-  - 


6. 1  Optimum  Thrust  Program 

In  Chapter  4  it  was  mentioned  that  if  the  maximum  flow  rale  of  a  rocket 
is  independent  of  position  then  the  adjoint  equations  become  uncoupled  from  the  equa¬ 
tions  of  motion  for  linearized  gravity  fields.  This  may  be  seen  by  writing  Eq.  (3.91 
for  this  case  as  Eq.  ((?.  I).  • 


(6.1) 


The  differential  equation  for  the  primer  vector  depends  only  upon  the  time  history  of 
the  gravity  gradient.  In  a  linearized  analysis  the  change  in  the  gravity  gradient  due 
to  a  change  in  position  is  neglected.  Therefore,  the  adjoint  equations  for  the  primer 
vector  become  uncoupled  from  the  equations  of  motion.  The  gravity  vector  and 
gravity  gradient  matrix  arc  evaluated  along  a  nominal  trajectory.  As  long  as  position 
deviations  from  the  nominal  are  small,  such  an  analysis  will  give  a  good  approximation 
to  the  motion. 

The  nominal  trajectory  may  be  powered  or  unpowered.  A  particularly 
important  case  is  where  the  nominal  trajectory  is  an  unpowered  elliptic  orbit.  Optimal 
low-thrust  trajectories  in  the  vicinity  of  such  an  elliptic  orbit  may  be  found  analytically 
by  lineai'izing  around  this  elliptic  orbit.  Since  the  primer  vector  history  along  a  coast¬ 
ing  arc  was  found  in  the  previous  chapter  this  primer  vector  history  will  constitute  the 

optimal  control  for  the  linearized  problem. 
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One  consequence  of  the  linearization  will  be  that  the  Hamiltonian  for  the 
power  limited  problem  will  no  longer  be  a  constant  of  the  motion.  The  reason  for  this 
is  that  the  gravity  gradient  matrix  in  Eq.  (G.  1)  is  now  a  function  of  time.  In  order  to 
have  a  constant  Hamiltonian  it  would  be  necessary  to  solve  the  exact  nonlinear  equa¬ 
tions  of  motion.  The  fact  that  the  Hamiltonian  is  not  a  constant  for  the  linear  problem 
creates  no  difficulty  as  a  full  set  of  integrals  for  this  problem  is  obtained  during  the 
course  of  the  explicit  solution  of  the  problem. 

The  linearized  motion  in  the  vicinity  of  an  elliptic  orbit  is  conveniently 
solved  in  a  variation  of  parameter  formulation.  The  equations  of  motion  are  written 
as  Eqs.  (0.2)  to  (6.  7). 


da 

dt 


=  2 


F  ef  +  Fpf 
r  1  p  5 

niL 
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dt  J  y  mfg 


(6.3) 


d0  = 
dF  = 
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dt 


(6.6) 


F  f 
2  1 


ill 


dM 

dt 


i>  -W6<f3  +  4» 


mcf 


(6.7) 


The  angular  variables  0,  a,  and  8  represent  small  rotations  of  the 
elliptic  orbit  in  its  own  plane,  around  the  major  axis,  and  around  the  latus' rectum, 
respectively.  As  the  linearization  of  the  problem  only  involves  small  rotations  it  is 
convenient  to  use  orthogonal  rotations  rather  than  Euler  angles.  The  angular  variables 
are  defined  by  Eqs.  (6.8)  through  (6.10)  in  terms  of  variations  in  the  conventional  Euler 
angles. 

0  s  8  oj  +  cos  i  6  0 

a  s  cos  W  6  i  +  sin  i  sin  o>  60 

8  =  -sin  w6  i  J-  sin  i  cos  co  6  O 


(G.  81 


(6.9) 


(6. 10) 


The  Lagrange  multipliers  for  the  linear  problem  will  be  exactly  the 
same  as  for  the  coasting  arc.  That  is,  the  Lagrange  multipliers  for  each  of  the 
orbit  elements  except  the  semimajor  axis  will  be  constant.  The  Lagrange  mul¬ 
tiplier  for  the  semimajor  axis  will  be  a  linear  function  of  time  given  by  Eq.  (6. 11) 
and  evaluated  on  the  unpowered  coasting  ellipse.  / 
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(G.  11) 


For  this  set  of  orbit  elements  the  optimal  thrust  acceleration  program  is  given  bj'  Eqs. 
(G.  12)  through  (G.  14). 


2Xaae  t,  +  Xeef/5  -  (XQ  -  fgX  M)  f4f5  ~  2XMC13 


ef. 


<6. 12) 


m  xe 


P  2Xaaef5  +  XeCVf4  *  W  +  <X0  “  f5V<f3  +£>f! 
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ef„ 


(G.  13) 


Jy  f5 


(G.  14) 


Thr  thrust  acceleration  program  represents  (he  vector  sum  of  six  different  accelera¬ 
tion  programs  corresponding  to  (he  initial  values  of  each  of  the  Lagrange  multipliers. 
As  the  thrust  acceleration  vector  for  this  problem  is  identical  with  the  primer  vector 
on  a  coasting  ellipse,  these  acceleration  programs  are  also  of  intei’est  for  impulsive 
transfers.  The  programs  corresponding  to  the  first  five  Lagrange  multipliers  are 
illustrated  in  Figures  (G.  1)  through  (0.5)  for  an  eccentricity  equal  to  the  square  root 
of  1/2.  Each  of  these  programs  will  maximize  the  change  in  its  orbital  element  for 
a  given  fuel  consumption  during  a  fixed  time  period.  In  general,  each  program  will 
also  produce  changes  in  the  other  orbital  elements. 
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The  optimum  program  for  changing  semimajor  axis  is  illustrated  in 
Figure  (0.1).  For  this  program  the  primer  vector  is  always  directed  along  the 
velocit}'  vector  and  its  magnitude  is  proportional  to  the  magnitude  of  the  velocity. 

The  locus  of  the  tip  of  the  primer  vector  will  describe  a  circle,  as  this  circle  is  known 
to  be  the  hodograph  of  the  velocity  vector  for  Keplerian  motion.  The  primer  vector  will 
have  a  stationary  maximum  at  periapsis  and  a  minimum  at  apoapsis.  At  an  eccentricity 
of  zero,  the  primer  vector  will  remain  constant  in  magnitude  as  it  moves  around  the 
circular  orbit.  As  the  eccentricity  goes  to  unity,  the  magnitude  of  the  primer  vector 
at  periapsis  will  approach  infinity.  If  this  magnitude  is  constrained  to  be  finite  then 
the  primer  vector  will  consist  of  a  discontinuous  function  having  finite  magnitude 
at  periapsis  and  being  zero  everywhere  else. 


Fig.  (G.  1)  Optimum  thrust  program  for  change  in  a(c  -  0,  707). 
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The  optimum  thrust  program  for  changing  eccentricity  is  illustrated 
in  Figure  (G.2).  This  program  always  has  two  equal  maxima  at  apoapsis  and  at  peri- 
apsis  for  all  eccentricities.  The  primer  vector  solution  corresponding  to  this  Lagrange 
multiplier  on  a  coasting  arc  allows  for  two  impulse  transfers  with  this  coasting  arc  in 
between  the  two  impulses.  The  direction  of  the  thrust  always  approximates  the  direc¬ 
tion  of  the  semiminor  axis  of  the  ellipse.  For  small  eccentricities,  there  is  1  latively 
little  variation  in  the  magnitude  of  the  acceleration.  At  an  eccentricity  of  unity,  the 
magnitude  of  the  acceleration  vector  varies  linearly  from  zero  at  the  origin  to  its 
maximum  value  at  apoapsis.  Like  the  semimajor  axis  program  at  unit  eccentricity 
there  is  also  an  isolated  value  at  periapsis  which  in  this  case  is  equal  to  the  value 
of  the  primer  vector  at  apoapsis.  For  this  eccentricity  of  unity,  the  thrust  is  direc¬ 
ted  exactly  at  right  angles  to  the  major  axis. 


Fig.  (6.2)  Optimum  thrust  program  for  change  in  e(e  -  0.707). 
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\  The  optimum  thrust  program  for  rotating  the  orbit  in  its  own  plane  is 

shown  in  Figure  (G.3).  In  this  ease,  the  primer  vector  also  has  two  equal  maxima 

which  constitute  possible  locations  for  impulses  in  the  impulsive  problem.  These 
maxima  occur  in  the  vicinity  of  the  semiminor  axis  and  arc  exactly  at  the  semiminor  axis 

for  eccentricities  of  zero  and  one.  For  an  eccentricity  of  zero,  this  program  is  the  same 

as  the  program  for  changing  eccentricity  but.  is  moved  90°  around  the  orbit.  For  an 

eccentricity  of  unity  the  thrust  is  directed  at  right  angles  to  the  major  axis  instead  of 

lying  close  to  it.  At  this  eccentricity,  the  acceleration  program  varies  smoothly  from 

zero  at  apoapsis  and  periapsis  to  a  maximum  at  the  semiminor  axis.  At  this  eccentricity, 

this  program  is  the  same  ns  the  program  for  rotating  the  orbit  around  the  latus  rectum 

except  that  it  lies  in  the  plane  of  the  orbit. 


Fig  (6.3)  Optimum  thrust  program  for  change  in  9  (e  ~  0.707). 
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Fip.ro  (0.4)  shows  the  optimum  program  for  rotating  the  orbit  around  the 
semimajor  axis.  fa  this  case,  the  magnitude  of  the  primer  vector  is  proportions  to  its 
perpendicular  distance  from  the  latus  rectum.  Figure  (6.5)  shows  the  primer  vector 
solution  for  rotating  the  orbit  around  the  la, us  rectum.  this  case,  the  magnitude 

of  the  primer  vector  is  proportional  to  the  perpendicular  distance  from  the  maj  or 
axis. 


Tig.  (G.4)  Optimum  thrust  program  for  change  in  ol  (e  =  0.  707). 


Fig,  (0.5)  Optimum  thrust  program  for  change  in  /3  (c  =  0.707). 


The  sixth  thrust  program,  (he  program  for  changing  (he  mean  anomaly, 
is  a  combination  of  three  programs,  the  program  for  rotating  the  orbit  in  its  own  plane, 
a  purely  radial  program  whose  magnitude  is  proportional  to  radius  and  the  semimajor 
axis  program  with  a  linear  lime  variation  superimposed  upon  it.  Because  of  the  chano-e 

in  the  semimajor  axis  program  with  lime,  this  is  the  only  program  that  will  not  be  peri- 
odic  and  will  not  repeat  itself. 
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6.2  Integration  of  the  Equations 


The  correct  first  order  terms  for  the  rates  of  change  of  all  elements 
except  the  mean  anomaly  may  be  found  by  simply  substituting  Eqs.  (6. 12)  through 
(6.14)  into  Eqs.  (6.2)  through  (6.6)  and  integrating.  For  these  integrations  the  orbit 
elements  are  taken  as  the  elements  of  the  unperturbed  ellipse.  The  independent  vari- 

'  t 

able  used  for  the  integrations  is  the  eccentric  anomaly  of  the  unperturbed  orbit.  The 
first  order  perturbations  in  semimajor  axis  produce  first  order  perturbations  in  mean 
anomaly  so  lhatEq.  (6.15)  must  be  used  to  calculate  the  perturbations  in  the  mean 
anomaly.  The  double  integration  implied  by  this  equation  may  be  avoided  by  rewrit¬ 
ing  the  perturbation  in  semimajor  axis  in  terms  of  the  identity  of  Eq.  (6. 16). 


dM  y  3 

F.  F* 

dt  J  a3  2, 

I a  +  / 

/a  Jy  mef3 

_  dltinj  _  da 
dt  dt 

f6, 15) 


(6.16) 


The  integration  of  the  equations  of  motion  with  the  optimum  thrust  pro¬ 
gram  is  straight  forward  and  results  in  Eqs.  (6, 17)  through  (6.24).  The  various  inte¬ 
grals  occuring  in  these  integrations  arc  expressed  by  the  <£  functions  given  following 
Eq.  (6.24). 
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These  equations  represent  a  complete  first  order  solution  for  optimum 
transfers  in  the  neighborhood  of  an  elliptic  orbit.  They  are  linear  in  both  the  changes 
and  in  the  slate  and  the  corresponding  Lagrange  multipliers.  By  using  the  proper 
boundary  conditions  or  transversality  conditions,  it  is  possible  to  solve  any  desired 
interception, rendezvous  or  orbit  transfer  problem.  For  interception  the  transvers¬ 
ality  condition  is  that  the  final  value  of  the  thrust  acceleration  given  by  Eqs.  (6.12) 
through  (6.14)  must  be  zero.  For  rendezvous  the  initial  values  of  all  the  Lagrange 
multipliers  must  be  found  so  as  to  drive  the  terminal  state  deviations  to  zero.  For 
orbit  transfer,  the  transversality  condition  is  that  the  Lagrange  multiplier  for  the 
mean  anomaly  must  be  zero.  The  other  five  Lagrange  multipliers  must  be  deter¬ 
mined  so  as  to  reduce  the  deviations  in  their  five  state  variables  to  zero. 

6. 3  Secular  Changes  in  the  Orbit 

While  Eqs.  (6.17)  through  (6.24)  represent  a  complete  solution  of  the 
problem  they  are  rather  complex  and  contain  a  large  number  of  terms.  If  the  time 
for  rendezvous  or  orbit  transfer  is  large,  so  that  the  maneuver  requires  many  rev¬ 
olutions  in  the  elliptic  orbit,  then  the  equations  can  be  greatly  simplified  and  an  ex¬ 
plicit  solution  is  easily  obtained.  This  is  done  by  neglecting  the  bounded  periodic 
terms  in  Eqs.  (6. 17)  to  (6.24)  in  comparison  with  the  dominant  secular  terms  to 
yield  Eqs.  (6.25)  to  (6.32). 
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For  rendezvous  or  orbit  transfer,  the  Lagrange  multipliers  that  satisfy  the  boundary 
value  problem  are  easily  determined  and  the  payoff  may  be  written  directly  in  terms 
of  the  changes  in  the  elements  and  the  total  transfer  time. 
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If  only  long  time-orbit  transfer  is  considered,  all  the  cross  product 
terms  in  Eqs.  ((>.25)  through  (6. 32)  disappear.  The  thrust  program  corresponding 
to  the  Lagrange  multiplier  for  each  orbit  element  not  only  maximizes  the  rate  of  change 
of  that  element,  but  produces  no  change  in  the  other  four  elements.  This  particular  set 
of  orbit  elements  was  chosen  because  they  have  this  orthogonality  property.  Not  all 
sets  of  orbit  elements  will  have  this  property. 

Equations  (6. 12)  through  (6. 141  show  that  there  are  no  secular  changes 
in  the  thrust  acceleration  programs  for  orbit  transfer.  The  averaged  thrust  accelera¬ 
tion  will,  therefore,  remain  constant  from  revolution  to  revolution,  and  the  orbit  ele¬ 
ments  will, on  the  average,  change  linearly  with  time. 

For  rendezvous  problems  the  final  value  of  the  mean  anomaly,  the  sixth 
orbit  element  must  be  considered.  Equation  (6.33)  represents  the  fuel  requirement 
for  changing  all  six  elements  over  long  time  periods.  The  quantity  Am  defined  by 
Eq.  (6.34)  represents  the  perturbation  in  the  mean  anomaly  produced  by  the  thrust 
programs  for  changing  semimajor  axis  and  for  rotating  the  semimajor  axis  in  its 
own  plane.  The  fuel  required  to  change  the  mean  anomaly  depends  upon  the  square 
of  the  difference  between  the  desired  perturbation  and  this  particular  perturbation. 

The  thrust  program  that  produces  the  desired  perturbation  in  the  mean  anomaly  is 
the  semimajor  axis  program  of  Figure  (6.1)  with  an  averaged  thrust  acceleration 
that  is  a  linear  function  of  lime.  The  tenns  involving  \  in  Eq.  (6. 12)  and  (G.  13) 
have  a  neglignble  effect  on  long-time  motion  but  are  significant  for  short  time  ren¬ 
dezvous.  Since  Eq.  (6.33)  holds  only  in  the  case  of  long  time  motion,  the  mean 
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anomaly  term  of  Eq.  (6.33)  will  generally  be  negligible  compared  to  the  terms  in 
bracket?  because  it  is  inversely  proportional  to  time  cubed.  The  fuel  required  for 
rendezvous  will  be  only  slightly  greater  than  the  fuel  reqxiircd  for  orbit  transfer  for 
these  long-time  cases. 
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CHAPTER  7:  OPTIMUM  POWER-LIMITED 
ORBIT  TRANSFER  IN  STRONG  GRAVITY  FIELDS 

7.1  Integrals  of  Motion 

In  section  6.3  of  the  previous  chapter,  the  secular  rates  of  change  of  the 
orbit  elements  in  the  presence  of  small  perturbing  thrusts  were  determined.  In  the  pre¬ 
sent  chapter,  these  rates  of  change  of  tine  orbit  elements  will  be  integrated  over  many 
revolutions  to  determine  the  optimum  sequence  of  orbit  elements  between  the  initial 
orbit  and  the  final  orbit.  The  changes  in  the  orbit  elements  will  no  longer  be  assumed 
to  be  small,  but  the  rates  of  change  of  the  orbit  elements  due  to  a  small  perturbing 
thrust'  will  be  assumed  to  be  small.  Only  the  first  order  terms  in  the  equations  of 
motion  will  be  considered.  The  analysis  is  an  application  of  what  is  known  as  Kryloff 
Bogoliuboff  averaging.  The  errors  of  this  type  of  averaging  will  be  on  the  order  of 

the  square  of  the  ratio  of  the  thrust  acceleration  to  the  acceleration  of  gravity.  For 

-4 

typical  electric  propulsion  systems  with  accelerations  of  about  10  g's,  the  error 
in  this  approximation  will  be  on  the  order  of  1  percem  at  about  10  earth  radii. 

Using  the  l’esults  of  section  G.  3  of  the  previous  chapter  an  averaged 
Hamiltonian  may  be  written  as  Eq.  (7. 1). 


II  =■ 
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8  2  6  2e2  9 


+  -Ltis  >241x2 
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(7.1) 


For  this  Hamiltonian  the  Lagrange  multiplier  for  the  mean  anomaly  has  been  assumed 
to  be  zero  ns  only  orbit  transfer  and  not  rendezvous  will  be  treated.  The  Lagrange 
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multipliers  for  the  angular  variables  in  Eqs.  (7. 1)  may  be  expressed  in  terms  of  the 
Lagrange  multipliers  of  the  conventional  Euler  angles  by  Eqs.  (7.2),  (7.3),  and  (7.4b 

^0  ^  W  (7.2) 

=  fii  ~  X.  cos  ton-  X^sin  to  esc  i  -  X  sin  OJ  cot  i  (7.3) 


X r>  -  &n  -  -X.  sin  UH-  X  „  cos  to  esc  i  -  X  cos  to  cot  i  (7.4) 

P  2  l  Q  to 

I 

The  use  of  the  conventional  Euler  angles  will  allow  large  rotations  to  be  treated.  It 

should  be  noted  that  X  and  \Q  are  the  same  as  the  quantities  &  and  8  of  Eq.  (5.32) 

a.  P  1  2 

The  vector  integral  given  by  Eqs.  (5.44)  through  (5.48)  of  Chapter  5  apply  also  for  the 
present  problem. 

It  is  convenient  to  replace  the  eccentricity  by  a  new  variable  o  defined 
as  the  arc  sign  of  the  eccentricity,  Eq.  (7.5).  Its  corresponding  Lagrange  multiplier 
is  given  by  Eq.  (7.  G) 


In  terms  of  this  new  variable,  the  averaged  Hamiltonian  is  now  given  by  Eq.  (7.7). 
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(7.7) 
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It  should  bo  noted  that  the  l'ates  of  change  of  the  semimajor  axis  and  the  variables  <p 
and  0  are  independent  of  the  eccentricity.  Only  the  rates  of  change  of  the  rotations 
around  the  major  axis  and  in  the  orbit  plane  depend  upon  the  eccentricity.  By  averag¬ 
ing  the  Hamiltonian  of  Chapter  5,  another  integral  given  by  Eq.  (7.  8)  is  obtained. 


II 


=  const. 


(7.8) 


This  integral  may  also  be  obtained  by  noting  that  the  Hamiltonian  is  independent  of 
the  independent  variable,  time.  This  implies  that  the  averaged  value  of  the  square 
of  the  thrust  acceleration  is  a  constant  throughout  the  motion.  Equation  (7. 8)  may  be 
immediately  integrated  to  determine  an  integral  for  the  payoff  given  by  Eq.  (7.9). 


J  =  Ht  (7-9) 

Because  the  averaged  acceleration  is  a  constant,  the  cost  increases  linearly  with  time. 
Another  integral  may  be  found  by  averaging  Eq.  (5. 18)  and  (5.49)  to  yield  the  results 
given  by  Eq.  (7.30). 
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Another  integral  may  be  found  by  writing  down  the  rates  of  change  of  the 
semimajor  axis  and  its  Lagrange  multiplier.  The  averaged  Hamiltonian  possesses  the  con¬ 
ventional  canonical  equations  in  terms  of  the  averaged  values  of  the  orbit  elements  and 
their  Lagrange  multipliers. 


£-*x  V 

dt  a  y  (7.ii) 
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By  dividing  Eq.  (7.12)  by  (7.11)  these  equations  may  be  integrated  to  yield  Eqs.  (7. 13) 
and  (7.14) 
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Equation  (7. 10)  may  be  verified  directly  from  these  equations.  The  notation  x°n  is 
defined  as  the  initial  value  of  x  raised  to  the  n  th  power. 


The  explicit  solution  for  the  semimajor  axis  and  its  Lagrange  multi¬ 
plier  allows  the  elimination  of  these  variables  and  the  reduction  of  the  problem  from 
a  five-dimensional  problem  to  a  four-dimensional  problem.  In  order  to  do  this,  it 
is  convenient  to  rewrite  Eq,  (7. 7)  in  terms  of  the  rates  of  change  of  the  orbital  ele¬ 
ments  rather  than  in  terms  of  the  Lagrange  multipliers.  The  equations  of  motion 
such  as  Eq.  (7. 11),  may  be  solved  for  their  Lagrange  multipliers  in  terms  of  the 
rate  of  change  of  the  corresponding  orbital  element.  This  can  easily  be  done  in 
the  present  case  because  of  the  lack  of  cross-product  terms  to  yield  Eq.  (7. 15). 


2a  4a 


*  2 

a  2?2 
—  +  -*  + 


1  +  5  cot 


1  +  5  tan 


(7.15) 


A  new  payoff  variable  if)  whose  rate  of  change  is  defined  by  Eq,  (7. 16)  is  now 
introduced. 
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1  +  5  cot 
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(7.16) 


It  will  now  be  possible  to  separate  the  original  optimization  problem  into  two  parts. 
The  first  part  is  to  determine  the  cost  for  given  changes  in  the  semimajor  axis  and 
i}>  .  The  second  part  will  be  to  minimize  i|>  for  given  changes  in  the  remaining  four 
orbit  elements.  Eq.  (7. 16)  shows  that  this  second  problem  may  be  interpreted  geo¬ 
metrically  as  the  problem  of  determining  a  minimum  length  trajectory  or  geodesic 
in  a  four-dimensional  orbit  element  space. 
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In  terms  of  this  new  variable,  the  Hamiltonian  of  Eq.  (7. 7)  is  now  given 

by  Eq.  (7.17). 
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The  rate  of  change  of  ip  is  given  by  Eq.  (7. 18)  where  it  should  be  noted  that  \ip  is  a 
constant. 
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Equation  (7. 18)  may  be  integrated  by  use  of  Eq,  (7. 12)  and  (7. 13)  to  yield  Eq.  (7. 19). 
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The  constants  of  Eqs,  (7. 19)  and  (7. 14)  may  be  eliminated  to  determine  the  payoff 
explicitly  in  terms  of  the  initial  and  final  values  of  the  semimajor  axis  and  the  total 
change  in  ip  ,  Eq.  (7.20). 
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This  equation  has  an  interesting  physical  interpretation.  The  cost  is  the  same  as  the 
cost  of  transfering  in  field  free  space  from  the  mean  orbital  velocity  of  the  initial  orbit 
to  the  mean  orbital  velocity  of  the  second  orbit  with  an  angle  of  the  between  the 
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two  velocity  vectors.  By  the  definition  of  0  this  will  provide  a  complete  solution  for 
problems  where  the  eccentricity  is  changed  but  there  is  no  rotation  of  the  elliptic  orbit. 
In  such  problems,  the  variable  0  will  be  equal  to  the  total  change  in  the  orbit  element 
<p  divided  by  i/P. 


105 


7.2  Coplanar  and  Coaxial  Transfers 


It  is  possible  to  solve  the  four-dimensional  problem  of  minimizing  0  in 
terms  of  the  remaining  four-orbital  elements  explicitly.  However,  the  solution  of  this 
general,  orbit  transfer  problem  is  rather  complex  and  only  2  two-dimensional  problems 
will  be  considered  in  this  chapter.  These  two  problems  ax'C  the  important  eases  of  co- 
planar  and  coaxial  transfers.  In  both  cases  there  is  only  a  rotation  around  a  single  axis 
In  the  coplanar  case,  the  rotation  is  around  an  axis  perpendicular  to  the  orbit  plane.  In 
the  coaxial  case  the  rotation  is  around  the  major  axis  of  the  orbit.  The  coplanar  and 
coaxial  cases  arc  obtained  from  the  general  case  by  means  of  the  following  definitions. 

Coplanar 

0  =  to 

i  =  0  xc,a° 

0=0  ~  0 

Coaxial 


to  =  0 


i  =  a 


0=  0 


0 
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Because  of  the  particular  form  of  the  Hamiltonian,  Eq,  (7.7),  the  coplanar  and  coaxial 

v 

cases  can  be  transformed  into  one  another  so  that  it  is  only  necessary  to  solve  one  of 
these  problems  to  obtain  the  solution  of  both.  The  following  transformation  equations 
will  take  the  solution  of  one  problem  into  the  solution  of  the  other  problem. 
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Coaxial 
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(7.28) 


The  variable  k  in  Eq.  (7.28)  is  defined  by  Eq.  (7.29). 


n2  .2  ■.  .  02  „  2  2  0 

A  cot  k  “  X  +  X  col  o 

03  1  <P  CU 


(7.29) 


By  utilizing  the  definition  of  the  rale  of  change  of  &  given  in  Eq.  (7. 16)  and  integrating. 
Eq.  (7.30)  for  0  is  obtained. 


(7.30) 

'  /• 
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Equations  (7.28) 
anti  to  in  terms 


and  (7.20)  provide  a  parametric  representation  of  the  solution  for 
of  i!)  and  the  initial  conditions  on  p  and  the  Lagrange  multipliers. 


$ 


0 

o  ) 


(7.31) 


p(4>, 


(7.32) 


This  type  of  representation  can  be  explicitly  solved  for  in  the  present  case.  The  struc¬ 
ture  of  the  resulting  extremal;  is  shown  in  Fig.  7. 1.  The  values  of  u>  are  limited  to 
the  range  between  0  and  tt.  Larger  rotations  need  not  be  considered.  The  values  of 
p  are  limited  to  the  range  between  0  and  2TT.  Values  of  p-  greater  than  rr  represent 
elliptic  orbits  having  the  opposite  sense  of  rotation  to  the  initial  orbit.  The  direction  of 
rolnlion  is  reversed  by  passing  through  a  unit  eccentricity  ellipse  having  zero  angular 
momentum  but  a  finite  semimajor  axis. 


Fig.  7. 1 

Structure  of  the  extremals  for  e>0  -  45°. 


109 


If  the  desired  terminal  value  of  <p  is  also  00°  there  will  be  no  change  in  the  eccentri¬ 
city  during  the  transfer  if  the  total  change  in  the  angle  CO  is  less  than  30°.  The  verti¬ 
cal  line  for  k  ~  90°  in  Fig.  7.2  is  a  minimizing  extremal  until  it  becomes  tangent  to 
the  envelope  at  this  point.  Beyond  30°  the  ellipses  decrease  in  eccentricity  and  then 
incroase'it  to  produce  the  extremals  shown  for  other  values  of  k^.  While  this  trans¬ 
fer  is  somewhat  academic  the  corresponding  coaxial  transfer  is  the  more  interesting 
case  of  transfer  between  inclined  circular  orbits.  For  inclinations  of  less  than  36° 
the  orbital  transfer  will  involve  no  change  in  the  eccentricity  of  the  intermediate 
orbits  while  for  angles  greater  than  30° the  intermediate  orbits  will  become  eccentric. 
The  extremals  and  payoff  curves  for  two  other  initial  values  of  o  ,  45° 
and  0°  are  shown  in  Figs,  7.3  and  7,4. 


Fig.  7.3  Extremals  and  payoff  curves  for  coplanar  transfer  with  o0  -  40°. 


^.DEGREES 

Fig.  7.4  Extremals  and  payoff  curves  for  coplanar  transfer  in  the  degenerate  case  oO  ---  0°. 
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In  Figure  7. 1  throe  sample  extremals,  A,  B,  and  C,  arc 'shown  representing  different 
values  of  the  constant  k  for  a  value  of  <p°  of  45°.  If  extremal  A  is  followed  from 
its  start  at  0  =  0,  it  will  first  pass  through  a  region  in  which  there  are  no  other  ex¬ 
tremals  which  go  to  the  same  point.  Up  to  point  1,  extremal  A  will  represent  a  unique 
solution  of  the  boundary  value  problem  for  any  point  it  passes  through.  At  point  1,  ex¬ 
tremal  A  first  crosses  an  envelope  of  extremals  shown  as  the  dotted  line.  Beyond  point 
1,  there  is  another  family  of  extremals  such  as  extremal  C  which  goes  to  every  point 
along  extremal  A.  The  boundary  value  problem  no  longer  has  a  unique  solution.  Be¬ 
tween  point  1  and  point  2,  extremal  A  represents  the  absolute  minimum  value  of  $ 
to  reach  any  point  along  it.  At  point  2,  extremal  A  has  the  same  cost  as  extremal 
C.  A  point  where  two  extremals  of  different  families  have  the  same  cost  is  known  as  a 
Darboux  point.  Between  points  2  and  3,  extremal  A  still  represents  a  locally  mini¬ 
mizing  solution  of  the  optimization  problem  but  does  not  represent  the  absolute  mini¬ 
mum  which  occurs  with  extremals  of  the  other  family.  At  point  3,  extremal  A  first 
becomes  tangent  to  the  envelope  of  extremals.  Such  a  point  is  called  a  conjugate  point. 
At  the  conjugate  point,  the  Jacobian  of  the  state  with  respect  to  the  initial  values  of  the 
Lagrange  multipliers  becomes  equal  to  0. 


t&l ija  « 
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(7.33) 
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If  extremal  A  is  continued  beyond  point  3,  it  ceases  to  be  even  locally  minimizing. 

That  is,  there  arc  curves  in  the  immediate  neighborhood  of  extremal  A  which  will 
have  a  lower  cost.  At  point  3,  extremal  A  violates  an  additional  neccssni'y  condi¬ 
tion  for  optimality  of  the  classical  calculus  of  variations  that  is  known  as  the  Jacobi 
condition.  The  Jacobi  condition  is  simply  that  a  minimizing  extremal  must  possess 
no  conjugate  points  between  its  initial  and  final  points.  For  problems  in  the  classi¬ 
cal  calculus  of  variations  where  the  control  never  lies  on  a  boundary,  the  combination 
of  the  strong  forms  (without  equality)  of  the  maximum  principle  and  the  Jacobian  condition 
is  sufficient  for  an  arc  to  be  locally  minimizing.  It  should  be  noted  that  these  suffici¬ 
ency  conditions  do  not  guarantee  that  the  extremals  represent  an  absolute  minimum 
as  is  shown  by  the  example  of  extremal  A  between  points  2  and  3.  On  extremal  C, 
the  three  points  1,  2,  and  3  become  coincident. 

A  set.  of  extremals  and  constant  cost  lines  for  the  minimizing  extremals 
is  shown  in  Fig.  7.2  for  an  initial  value  of  o  of  90°. 


Fig.  7.2  Extremals  and  payoff  curves  for  coplanar  transfer  with  O0  »  90°. 


112 


7-*’  "’°  >m,Wcra  bocomos  degenerate  because  the  argument  „r  perigee  of  » 
circular  orlm  is  undefined.  In  this  ease,  the  fuel  consumption  depends  only  upon 
the  changes  in  eccentricity  and  semimajor  nsts  between  the  initial  and  final  orbits. 
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CHAPTER  8:  LINEARIZED  IMPULSIVE  TRANSFER 
IN  THE  VICINITY  OF  A  CIRCULAR  ORBIT 

8, 1  The  Primer  Vector  on  a  Circular  Orbit 

The  conventional  orbit  elements  used  in  Chapter  15  to  derive  the  equations 
for  the  primer  vector  have  singularities  at  zero  eccentricity  and  zero  inclination.  It  is 
desirable  to  use  another  set  of  orbit  elements  which  is  well  behaved  for  this  case  to  rep¬ 
resent  the  primer  vector  on  a  circular  orbit.  One  such  set  of  orbit  elements  is  defined 
by  Eqs.  (  8. 1) 

f 

c  ==  e  cos  (w  +  fi)  e  =  e  sin  (W  +  £i) 

x  y 

i  i  5  i  cos  i  =  i  sin  D 

x  y 

M  s  M  +  W  +  ^  E  =  E  +  0)  +  H 

x  x 

For  the  particular  case  of  zero  eccentricity  the  rates  of  change  of  these 
oi'bit  elements  arc  given  by  Eqs.  (8, 2)  to  (8.  7). 


de 

dt 


F  ' 

—  2  cos  E 
m  x( 


(8.3) 


114 
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(8.4) 


The  primer  vector  in  terms  of  this  set  of  orbital  elements  is  given  by 
Eqs.  (8.8)  through  (8.  10)  and  the  rate  of  change  of  the  primer  vector  is  given  by  Eqs 
(8. 11)  through  (8. 13). 
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-(-X  sin  K  +  X  cos  E  + 
a  cx  x  eg  x 


(8.12) 


•  1 

X  ”  —(-X.  sin  E  -i-  X.  cos  E  ) 
/  a  i  xi  x 

x  y 


(8.13) 


As  for  the  elliptic  oi'bit,  all  the  Lagrange  multiplier's  will  be  constant  except  the  Lagrange 

I 

multiplier  for  the  semi  major  axis  which  will  vary  lincarily  with  time  or  with  the  central 
angle.  Eqs.  (8.14)  and  (8.15). 


x 


(8.141 


(8.15) 


It  is  somewhat  easier  to  visualize  the  behavior  of  the  primer  vector  by  transforming  the 
origin  of  the  polar  coordinate  system  to  yield  Eqs.  (8, 1G)  to  (8. 10).  - 
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Tins  form  for  the  primer  vector  on  a  circular  orbit  may  be  used  to  treat  transfers  be¬ 
tween  neighboring  near-circular  orbits.  The  geometric  properties  of  the  primer  vector 
are  of  considerable  importance  in  this  analysis  and  will  be  illustrated  in  the  followin'*- 


sections. 
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8.2  Timc-Opon  Case— Nodal  Solutions 

For  tho  time-open  case,  the  Lagrange  multiplier  for  the  mean  anomaly  is 
0  and  the  components  of  the  primer  vector  may  be  written  as  Eqs.  (8.20)  through  (8.22). 
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These  equations  represent  the  equations  of  an  ellipse  in  three  space.  This  ellipse  is 


formed  by  the  intersection  of  a  2:1  elliptical  cylinder  parallel  to  the  Xr  axis  and  a 
plane  which  passes  through  the  intersection  of  the  cylinder  axis  with  the  X^X^  plane. 

A  typical  case  is  illustrated  in  Fig.  8.  1  which  also  shows  the  projection  of  the  ellipse 
on  the  X^Xq  plane.  v 


Primer  Locus  Diagram 


I-ig.  8.1 
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For  transfers  with  more  then  one  impulse,  this  elliptical  primer  locus 
must  huve  more  lh;m  one  equal  maxima.  There  are  only  three  configurations  of  this 
elliptical  primer  locus  which  allow  the  primer  vector  to  have  more  than  one  maximum. 
The  first  configuration,  to  be  treated  in  this  section,  represents  a  family  of  solutions 
where  the  center  of  the  ellipse  is  located  at  the  origin.  (Fig.  8.2)  In  this  ease  the  two 
equal  maxima  occur  on  the  major  axis  of  ihc  ellipse  and  arc  separated  by  180  .  for  this 
case  the  components  of  the  primer  vector  at  the  second  impulse  will  be  equal  in  magni¬ 
tude  but  opposite  in  sign  to  the  components  of  the  primer  vector  at  the  first  impulse. 


M 

Primer  Locus  Diagram 
Fig.  8.2 

As  the  two  impulses  arc  separated  by  180°,  both  impulses  must  occur  on 
the  line  of  intersection  or  line  of  nodes  between  the  initial  and  final  orbit  planes.  If  tho\ 
do  not,  then  the  second  impulse  could  not  remove  all  of  the  inclination  between  the  two 
planes.  As  a  result,  this  ease  will  be  referred  to  as  Ihc  nodal  case.  The  total  &Y  ior 
the  nodal  case  is  given  by  Kq.  (8.23). 
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11  should  be  noted  that 'the  fuel  consumption  is  independent  of  (he  change  in  semimajor 
axis  for  this  case.  In  order  for  the  primer  vector  to  h  ve  two  equal  maxima  for  this 
case,  the  two  inequalities  given  by  Eqs.  (8.24)  and  ( 8 . 25)  must  be  satisfied.  ( 


A  a 

o 

a 


(8. 24' 


(8.25) 


For  the  equations,  the  x  axis  has  been  assumed  to  be  aligned  with  the  line  of  nodes. 

8.3  Time-Open  Case:  Nondegencratc  Solutions 

A  second  configuration  which  allows  two  equal  maxima  of  the  primer  vector 
corresponds  to  cases  where  the  ellipse  passes  through  the  Xp  rxis  and  the  primer  vectors 
again  lie  in  single  plane.  (Fig.  8.3) 
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Primer  Locus  Diagram 
Fig.  8.3 


This  case  will  be  referred  to  as  the  nondegenerate  case  as  the  fuel  consumption  will 
depend  upon  the  changes  in  all  of  the  orbital  elements.  For  this  case  the  radial  and 
out  of  plane  components  of  the  primer  vector  have  equal  and  opposite  values  at  the 
two  impulses,  while  the  circumferential  components  have  equal  values.  For  this 
case,  the  cnctrnl  angle  is  not  restricted  to  180°.  The  total  fuel  consumption  for 
the  nondegenerate  case  is  given  by  Eq.  (8. 20). 
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8.1  Time -Open.  Case;  Ri n;-;nl ;i v  Solutions 

The  third  configuration  of  (lie  primer  vector  locus  which  allows  for  more 
than  one  maxima  of  the  primer  vector  is  a  combination  of  the  two  previous  cases  where 
the  primer  locus  passes  through  the  Xg  axis  and  where  the  center  of  the  ellipse  is  lo¬ 
cated  at  the  origin  of  the  coordinate  system.  This  particular  locus  is  tilled  so  that 
the  primer  vector  has  unit  magnitude  at  all  points.  This  is  an  example  oT  a  singular 
solution  where  the  orbital  location  of  the  impulses  can  not  be  determined  from  the 
primer  vector  solution.  In  fact,  for  this  linear  problem,  the  singular  solution  is 
ununique  in  that  there  are  an  infinite  number  of  transfers  involving  different  numbers 
and  locations  of  impulses  all  of  which  have  the  same  total  fuel  consumption.  This 
particular  locus  (Fig.  8.4)  is  only  one  of  two  that  occur  for  circular  reference  orbits. 
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Primer  Locus  Diagram 
Fig.  8.4 


The  other  singular  locus  is  a  special  case  of  the  nondegenerate  ease  where  the  ellipse 
shrinks  to  a  point  located  at  plus  or  minus  one  on  the  Ag  axis.  For  the  singular  case 
illustrated  in  Fig.  8.4,  the  total  &V  is  given  by  Ecj.  (8.29). 


(8.29) 


Once  again,  there  are  two  inequalities  which  must  be  satisfied  to  distinguish  the  singular 
case  from  the  two  other  cases  which  can  occur. 

The  fuel  consumption  of  the  different,  cases  is  illustrated  in  Figs.  S.5 


through  8.10.  In  each  case,  the  total  change  in  semimajor  axis  and  eccentricity  are 
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" '  W"h  C<ml',U,  S  °r  vnlnc.  of  ll,o  ,„,a,  change  inclination.  A]|  ll|roc 

elements  ere  normal,.,  dividing  then,  ,,v  „,e  tote,  Av  assuming  ,  Br;lvHalional 

constant  or  unitv.  The  t„,„I  change  in  inclinnii,,,,  eon.onrs  rnngc  from  n  numerical 
vnlue  of  one  „t  the  origin  to  a  value  of  nor,,  at  the  edge.  of  the  square  figure..  Enel 
contour  represents  a  change  of  .05  in  the  normalise.!  inclination  change. 
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The  angular  variable  shown  to  the  right  of  each  figure  represents  the  angle  between  the 
line  of  nodes  and  the  desired  direction  of  the  change  in  eccentricity.  In  Fig.  8.5,  we 
have  the  coaxial  ease  in  which  the  change  in  eccentricity  is  alligned  with  the  line  of 
nodes  of  the  two  orbits.  In  this  case,  the  singular  ease  docs  not  arise  and  only  the 
nodal  case  and  the  nondegenerate  case  occur.  In  Fig.  8.  G,  thor  is  a  5°  angle  between 
the  eccentricity  change  vector  and  (he  line  of  nodes  and  the  singular  case  starts  to  ap¬ 
pear  at  the  left  and  right  sides  of  the  figure.  The  curves  for  both  the  singular  ease 
and  the  nodal  case  ai'e  indcpdndcnt  of  the  change  of  seroimajor  axis  and  are  straight 
lines  whereas  the  curves  for  the  nondegenerate  case  start  to  develop  some  curvature. 

As  the  angle  between  the  eccentricity  change  vector  and  the  line  of  nodes 
increases,  the  nodal  region  decreases  and  the  singular  region  grows  until  in  Fig.  8.10 
the  nodal  region  has  completely  disappeared. 
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8.5  'l  ime  Fixed  (’uses 


For  these  e:ises  we  must  return  to  Eqs.  (8.  1G)  through  (8.  18)  including 
the  Lagrange  multiplier  for  the  mean  anomaly.  When  this  term  is  included  the  elliptic 
locus  becomes  stretched  out  to  form  a  cycloid-dike  curve  in  three  dimensions.  It  is  as 
though  one  took  an  elliptic  coil  of  wire  and  stretched  it  out  uniformly  in  the  direction. 
In  this  case  it  is  possible  to  have  transfers  with  up  to  six  impulses  in  the  three  dimen¬ 
sional  case  and  up  to  four  impulses  in  the  two  dimensional  case.  There  is  a  general 
result  that  linear  transfers  require  a  number  of  impulses  no  greater  than  the  number 
of  states.  In  the  time-open  ease  this  means  that  transfers  between  neighnoring  near- 
circular  orbits  may  require  up  to  five  impulses.  However,  these  eases,  in  fact,  re¬ 
quire  no  more  than  two  impulses  even  in  the  singular  ease.  However  for  time-fixed 
rendezvous,  it  is  easy  to  construct  two  dimensional  transfers  which  use  the  maximum 
number  of  four  impulses.  Three  dimensional  transfers  which  use  the  maximum  num¬ 
ber  of  six  impulses  h:  ve  been  determined  by  Marec.  and  by  Breakwell. 

For  the  two  dimensional  ease,  Trussing  has  considered  the  four  impulse 
transfers  and  has  shown  that,  all  of  these  transfers  have  symmetry  properties.  Me  con¬ 
structs  these  transfers  by  taking  the  two  dimensional  cycloid-like  locus  and  drawing  a 
circle  which  is  tangent  to  the  locus  at  two  points.  As  long  as  these  tnngenoics  produce 
local  maxima  of  the  primer  vector,  he  has  determined  the  stationary  interior  maxima 
of  the  primer  vector.  The  next  intersections  of  the  circle  with  the  primer  locus  yield 
the  nonslationary  terminal  maxima  of  the  primer  vector.  Trussing  has  shown  that  these 
solutions  exist  with  both  odd  and  even  numbers  of  loops  of  cycloid-like  curves  included  jn 
the  circle.  The  single  loop  transfers  even  exist  for  transfer  angles  slightly  below  180°, 
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trussing  further  ordered  the  fixed  lime  rendezvous  between  ncMghboriug  circular 
orbKs.  He  found  that  these  rendezvouses  may  require  four  impulses,  three  impulses 
with  initial  of  final  coast  periods;  three  impulses,  two  impulses  with  an  initial  or  final 
const  period  and  two  impulses.  There  is  also  a  class  of  singular  cases  corresponding 
to  the  coplanar  singular  case  which  may  have  any  distribution  of  impulses  and  coasts. 
His  results  show  (hat  the  structure  of  the  minimum  fuel  rendezvous  solutions  in  even 
this  simple  case  is  very  complex.  The  extensions  of  tins  work  to  some  speoia!  three- 
dimensional  eases  by  Bvcakwell  and  by  Mnroc  show  that  the  extra  dimension  adds  pro- 
prolionally  to  the  complexity  of  the  solutions. 
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CHAPTER  9; 


OPTIMUM  IMPULSE  TRANSFER 


IN  AN  INVERSE  SQUARE  FIELD 


9.1  Time:  Open  Transfer  Between  Cophmnr,  Coaxial  El  lips  os 

There  are  four  different  classes  of  coplanar,  coaxial  transfer  problems  de¬ 
pending  upon  whether  the  ellipses  intersect  or  not  and  upon  whether  the  apoapsis  of 
the  ellipses  point  in  the  same  or  opposite  directions.  If  the  orbits  intersect,  then 
the  optimum  transfer  is  always  a  Hohmann  transfer  from  the  highest  apsis  to  the 
opposite  apsis.  For  the  intersecting  orbit  ease  one  impulse  will  always  be  an  ac¬ 
celerating  impulse  and  the  other  impulse  will  always  be  a  decelerating  impulse.  The 
primer  vector  solution  corresponding  to  this  transfer  is  the  solution  corresponding 
to  the  Lagrange  multiplier  for  the  eccentricity  in  Chapter  G.  These  two-impulse 
transfers  must  be  compared  with  a  transfer  that  goes  out  to  infinity  and  comes  back 
again.  If  the  axes  are  aligned,  then  the  optimum  transfer  via  infinity  will  be  via 
two  parabolas  which  arc  tangent  to  the  peri apses  of  the  two  ellipses.  Transfer  from 
the  first  to  the  second  parabola  is  via  an  infinitcssimal  impulse  at  infinity.  If  the 
orbits  are  not  aligned,  the  same  two  transfer  parabolas  to  infinity  arc  used  but  these 
parabolas  must  be  connected  by  an  elliptic  orbit  at  infinity.  Marchal  has  given  a 
number  of  results  which  indicate  whether  the  optimal  transfer  is  a  two  impulse  trans¬ 
fer  or  via  infinity. 

| 

If  the  orbits  arc  non-intersecting  and  the  axes  arc  aligned,  the  optimal  two- 
impulse  transfer  is  again  from  the  highest  apsis  to  the  opposite  apsis  of  the  second 
ellipse.  If  the  orbits  are  non-intersecting  and  the  apses  are  opposed,  then  the  Hohmann 
transfer  from  cither  apsis  of  the  outer  orbit  must  be  considered  to  determine  the  true 
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optimum  two-impulse  transfer.  The  primer  vector  solution  corresponding  to  the  non- 
intersecting  case  has  both  impulses  cither  accelerating  or  decelerating  the  velocity 
and  has  a  primer  vector  solution  corresponding  to  the  Lagrange  multiplier  for  the  semi- 
minor  axis.  Once  again  there  is  a  three  impulse  transfer  via  infinity  (the  bi-elliptic 
transfer)  if  the  axes  arc  aligned,  and  a  four  impulse  transfer  via  infinity  if  the  axes 
arc  not  aligned.  For  the  aligned  case,  the  optimal  transfer  is  always  through  infinity 
if  the  ratio  of  the  two  periapsis  radii  is  greater  than  11.938  and  is  always  two  impulse 
Ilohmann  if  the  ratio  of  these  periapsis  radii  is  less  than  nine. 

9.2  Time  Open  Transfer  between  Hon-Coplanar  Coaxial  Ellipses 

As  for  the  coplanar  case,  all  of  the  impulses  for  coaxial  transfer  should  be 
on  the  line  of  apsides  (which  is  also  the  line  of  nodes)  and  should  have  no  radial  com¬ 
ponents.  In  the  non-planar  case,  there  is  an  additional  type  of  transfer  possible  which 
does  not  occur  in  the  planar  case.  Namely,  a  three  impulse  transfer  where  the  third 
impulse  is  giving  at  a  finite  radius.  One  important  class  of  transfers  for  which  the 
finite  radius  three-impulse  transfer  is  always  optimal  is  for  transfer  between  circular 
orbits  of  the  same  radius  for  inclination  angles  between  zero  and  60.185  degrees. 
Marchal  has  once  again  given  a  number  of  rules  for  determining  whether  the  transfer 
should  be  a  two  impulse  transfer,  a  three  impulse  transfer  with  a  finite  periapsis  ra¬ 
dius  or  a  transfer  via  infinity.  For  all  of  the  transfer  most  of  the  inclination  change 
is  done  at  the  largest  apoapsis  radius.  The  other  impulses  always  lie  within  about 
six  degrees  of  the  orbit,  plane  of  the  initial  and  final,  ellipse.  The  optimum  inclination 
change  for  each  of  the  impulses  can  readily  be  determined  either  from  the  primer  vector 
solution  or  from  an  analysis  by  Sun  which  reduces  the  problem  to  the  well  known  spider 
fly  problem. 
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9.3  Time  Open  Tran:; fur  Between  Coplanar  Ellipses 

Time  open  transfer  between  coplanar  but  non-coaxial  ellipses  fall  into  three 
different,  categories.  Two  impulse  transfers,  three  impulse  transfers  with  all  im¬ 
pulses  occurring  at  finite  radii,  and  transfers  via  infinity.  The  three  impulse  trans¬ 
fers  occur  rarely  and  only  for  highly  eccentric  initial  and  final  ellipses  over  a  small 
range  of  orientation  angles.  The  two  impulse  transfers  arc  of  two  different  families. 
One  family  has  both  impulses  either  increasing  or  decreasing  the  energy.  For  this 
family  both  impulses  occur  on  one  side  ol'  the  major  axis  of  the  transfer  ellipse.  The 
limiting  members  of  this  family  are  the  Uohmann  transfer  whose  primer  vector  solu¬ 
tion  corresponds  to  the  Lagrange  multiplier  for  the  semi-minor  axis  and  which  lias  a 
transfer  angle  of  180°;  and  the  primer  vector  solution  corresponding  to  the  Lawdcn 

spiral  which  has  a  transfer  angle  of  0°.  The  other  family  of  two  impulse  transfers 

/ 

has  one  impulse  being  an  accelerating  impulse  and  the  other  impulse  being  a  deceler¬ 
ating  impulse.  For  this  family,  the  two  impulses  lie  on  the  opposite  sides  of  the 
semi-major  axis  of  the  transfer  ellipse  and  the  transfer  angle  is  always  in  the  vicinity 
of  180°.  The  limiting  members  of  this  family  of  two  impulse  transfers  arc  the  re¬ 
verse  Uohmann  transfer  corresponding  to  the  Lagrange  multiplier  for  eccentricity 
and  the  symmetric  transfer  corresponding  to  the  Lagrange  multiplier  for  changing 
the  argument  of  perigee. 


At  every  point  on  the  ellipse  there  is  only  a  small  range  of  angles  in  between 
the  local  horizontal  and  the  local  direction  of  the  velocity  vector  in  which  impulses 
may  be  directed  for  an  optimal  transfer.  The  limiting  values  of  this  set  are  dctcr- 
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mined  by  the  members  of  the  two  different  families  of  two  impulse  transfers  mentioned 
in  the  last  paragraph.  A  single  impulse  transfer  from  the  initial  ellipse  or  onto  the 
final  ellipse  I  rum  the  transfer  ellipse  must  lie  within  this  range  of  useful  angles.  For 
all  eccentricities  below  a  critical  eccentricity  of  0.i)25  there  is  a  positive  range  of 
useful  angles  at  every  point  on  the  ellipse.  Above  this  critical  eccentricity,  there  is 
a  range  of  positions  between  perinpsis  and  the  semi-minor  axis  for  which  there  is  no 
positive  range  of  useful  angles  and  which  cannot  lie  employed  in  optimal  transfer. 
Above  this  critical  eccentricity  it  is  possible  to  have  transfers  between  the  two  differ¬ 
ent  families  of  two  impulse  transfers  so  that  infinilessiinal  throe  impulse  transfers  in 
the  vicinity  of  these  highly  eccentric  ellipses  become  possible.  At  an  eccentricity  of 


unity,  the  range  of  forbidden  positions  where  no  optimal  impulse  may  be  applied  extend 

t. 

all  the  way  from  periapsis  to  the  semi-minor  axis. 


9.4  Time  Open  Transfer  Between  an  Ellipse  and  a  Hyperbola 

This  is  one  transfer  problem  which  is  completely  solved  in  the  time  open  case. 
Unfortunately  the  optimal  solution  requires  infinite  time.  This  optimal  solution  is 
useful  for  generating  the  optimal  finite  time  transfer  for  long  transfer  times  by  means 
of  a  perturbation  analysis.  The  optimum  time  open  transfer  requires  five  impulses. 
The  first  impulse  is  a  tangential  impulse  at  the  periapsis  of  the  initial  ellipse  which 
puts  the  vehicle  onto  an  escape  parabola.  The  second  impulse  is  an  infinitcssimal  im¬ 
pulse  at  infinity  on  the  escape  parabola  and  it  transfers  the  vehicle  onto  an  ellipse  at 
infinity  which  transfers  the  vehicle  to  a  different  position  on  the  celestial  sphere.  The 
third  impulse  is  an  infinitcssimal  impulse  at  infinity  and  it  transfers  the  vehicle  back 
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onto  :i  parabola  which  will  have  its  peri  apsis  radius  at  the  minimum  allowable  radius 
determined  by  the  planetary  radius  and  atmosphere.  At  the  periapsis  of  the  second 
pai.tbola,  a  second  finite  impulse  is  applied  tangentially  to  transfer  the  vehicle  onto 
an  escape  hyperbola  having  the  proper  energy  and  the  proper  direction  of  the  asymp¬ 
totic  velocity  vector.  The  fifth  and  final  impulse  is  an  infinilessimal  impulse  at  in¬ 
finity  on  the  escape  hyperbola  to  produce  the  correct  angular  momentum  vector  for 
this  hyperbola. 


9,5  Time  Open  Transfer  Between  Two  Hyperbolas 

The  optimum  time  open  transfer  between  two  hyperbolas  is  completely  solved 
only  in  the  case  where  the  radius  of  the  attracting  body  is  zero  so  that  impulses  may 
be  applied  at  the  center  of  attraction.  In  thin  ease  the  optimal  transfer  is  a  six  im¬ 
pulse  transfer  where  all  the  impulses  arc  of  infinilessimal  magnitude  and  the  trans¬ 
fer  time  is  infinite,  For  the  case  where  there  is  a  minimum  allowable  radius,  the 
transfer  problem  becomes  quite  complicated.  By  applying  infinilessimal  impulses 
at  infinity,  the  problem  is  reduced  to  the  problem  of  transfer  between  two  asymptotic 
velocity  vectors.  The  optimal  transfer  always  lies  in  the  plane  of  these  two  velocity 
vectors.  This  problem  has  been  partially  solved  by  Marchal  in  some  recent  papers. 
For  some  cases  the  optimal  transfer  requires  four  impulses  and  the  use  of  an  inter¬ 
mediate  ellipse  at  infinity.  This  type  of  transfer  requires  the  spending  of,  an  infinite 
amount  of  time  in  the  vicinity  of  the  planet.  The  other  optimal  transfers  only  require 
spending  a  finite  lime  in  the  vicinity  of  the  planet  and  involve  energies  which  arc  al- 
,  ways  above  the  parabolic  energy  level. 
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When  the  transfer  is  not  made  through  the  parabolic  level,  there  are  never 
more  than  two  finite  impulses.  If  there  is  more  than  one  finite  impulse,  at  least 
one  of  these  impulses  always  occurs  at  infinity.  The  other  impulse  will  generally 
occur  somewhat  above  the  minimum  allowable  radius.  It  will  only  occur  at  the  mini¬ 
mum  allowable  radius  if  a  transfer  between  the  two  asymptotic  velocity  vectors  may 
be  made  by  means  of  a  tangential  impulse  at  the  minimum  allowable  radius.  If  this 
is  not  the  case,  the  impulse  at  a  finite  radius  will  not  be  tangential.  The  transfer 
hyperbolas  may  or  may  not  descend  to  the  minimum  allowable  radius. 

9.6  Coplannr  Time-Open  Angle-Open  Transfer 

Tliis  is  one  transfer  which  has  been  completely  solved  analytically  for  the 
two  impulse  case.  The  problem  is  stated  as  optimal  two  impulse  transfer  between 
two  different  radii  and  between  fixed  radial  and  circumferential  components  of  ini¬ 
tial  and  final  velocity.  It  may  be  looked  upon  as  transfer  between  specified  locations 
on  initial  and  final  orbits  with  the  argument  of  periapsis  of  the  orbits  being  left  open. 
This  solution  is  useful  for  re-entry  and  for  ascent  trajectories  and  has  guidance  ap¬ 
plications  since  it  can  be  written  down  in  closed  form.  In  general  these  transfers  can 
be  improved  by  adding  an  extra  degree  of  freedom  and  not  specifying  the  points  of  ar¬ 
rival  and  departure  on  the  terminal  orbits.  In  this  latter  case,  the  optimal  orientation 
is  always  coaxial  with  the  periapsis  pointing  in  the  same  direction. 

9.7  Time  Open  Transfer  Between  Fixed  End  Points 

The  problem  of  time-open  two-impulse  transfer  between  fixed  positions  on 
fixed  initial  and  final  orbits  has  been  considered  by  a  number  of  authors.  The  solution 


has  boon  reduced  to  finding  the  roots  of  an  eighth  degree  polynomial.  It  has 'been 
shown  that  for  some  initial  and  final  positions  on  elliptic  orbits,  the  optimal  two 
impulse  transfer  may  involve  a  hyperbolic  intermediate  orbit.  In  this  case,  ns  well 
as  in  some  other  cases,  the  fuel  consumption  of  the  two  impulse  transfer  may  be  re¬ 
duced  by  going  to  a  transfer  through  infinity.  If  the  constraints  that  a  finite  impulse 
must  be  applied  at  the  specified  initial  and  final  positions  is  dropped,  then,  by  allow¬ 
ing  coasting  in  the  initial  and  final  oi’bits,  this  problem  may  be  reduced  to  the  problem 
of  optimal  transfer  between  specified  initial  and  final  orbits.  The  problem  with  speci¬ 
fied  positions  does  have  some  application  to  ascent  and  descent  trajectories  where 
there  may  be  a  minimum  allowable  radius. 

9.8  Time  Fixed  Transfer 

Very  little  has  been  done  on  the  problem  of  time  fixed  transfer  and  rendez¬ 
vous.  The  work  that  has  been  done  has  shown  that  the  problem  is  one  of  considerable 
complexity.  The  work  with  neighboring  orbits  has  shown  that  the  optimal  transfer  may 
involve  up  to  six  impulses.  The  work  with  singular  arcs  has  shown  that  in  the  time 
fixed  case  there  may  exist  optimal  coplanar  singular  arcs  so  that  singular  arcs  as 
well  as  a  large  number  of  impulses  must  be  considered  in  finding  the  minimizing  ex¬ 
tremal  for  these  transfers.  The  major  contribution  has  been  an  iterative  technique 
which  allows  the  determination  of  an  optimal  n-impulse  transfer  from  the  non- optimal 
two-impulse  transfer  between  the  initial  and  final  positions.  However,  this  method 
is  only  a  method  for  (hopefully)  finding  a  locally  minimizing  solution.  In  general  it 
will  not  be  known  if  there  arc  other  n-impulse  transfers  or  if  there  are  other  trans- 


ar  arcs.  This  problem  is  a  fruitful  area 


CHAPTER  10:  SINGULAR  ARCS 


10.1  Necessary  find  Sufficient  Condi  lions  for  Singular  Arcs. 

For  optimal  rocket  trajectories  with  constant  exhaust  velocity,  a 
phenomenon  known  as  a  singular  arc  may  occur,  A  singular  arc  occurs  when  the 
coefficient  of  the  thrust  acceleration  in  the  Hamiltonian  remains  identically  zero 
over  a  finite  thrust  are  for  a  finite  time.  In  these  cases,  the  optimum  magnitude 
of  the  thrust  acceleration  can  not  be  determined  from  the  Hamiltonian,  Imposing 
the  condition  that  the  coefficient  of  the  thrust  acceleration  in  the  Hamiltonian  re¬ 
mains  identically  zero  will  allow  the  derivation  of  a  magnitude  for  the  thrust  ac¬ 
celeration  at  each  point  on  a  singular  arc.  For  the  case  where  the  thrust  accel¬ 
eration  may  become  unbounded,  this  will  correspond  to  the  magnitude  of  the  primer 
vector  being  unity  along  such  an  arc.  Considerable  information  about  such  arcs  may 

I 

be  determined  by  taking  the  successive  derivatives  of  the  magnitude  of  the  primer 
vector,  all  of  which  will  be  zero.  If  the  optimal  thrust  magnitude  is  bounded  then 
the  magnitude  of  the  primer  vector  will  be  constant  on  a  singular  arc.  As  has  been 
pointed  out  earlier,  for  linear  problems  singular  arcs  correspond  to  nonunique 
solutions  with  undefined  acceleration  magnitudes.  However,  for  the  nonlinear  prob¬ 
lem  the  acceleration  magnitude  on  the  singular  arcs  become  well  defined.  Such  arcs 
may  be  candidates  for  portions  or  all  of  a  minimizing  extremal. 

The  literature  on  the  classical  calculus  of  variations  provides  very 
little  information  on  singular  arcs.  In  recent  years  the  discovery  that  singular  arcs 
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mny  occur  in  minimizing  solutions  luts  lend  1o  renewed  interest  in  this  problem.  This 
lcccnt  work  hns  culm  inn  led  in  necessary  and  suliicient  conditions  for  singular  arcs 
although  Ihe  problem  of  coupling  singular  extremals  to  nonsingular  extremals  is  still 
not  in  a  satisfactory  stale.  The  first  of  the  new  necessary  conditions  was  generalized 
or  transformation  of  the  Legendre- Clebsh  condition  of  the  calculus  of  variations.  The 
extension  was  carried  out  by  n  number  of  investigators  including  Kelley,  Tait,  Robbins, 
Gurley,  Goh,  and  Speyer.  An  additional  necessary  condition  for  singular  arcs  was  then 
developed  by  Jacobson.  Move  recently  necessary  and  sufficient  conditions  have  been 
obtained  by  Jacobson,  Speyer  and  Jacobson,  and  McDanell  and  Powers.  /  These  condi¬ 


tions  allow  the  testing  of  a  singular  are  for  local  optimality. 


10.2  Junction  Conditions 

Kelley,  Kopp  and  Moyer  have  shown  that  optimal  junctions  between 
singular  and  nonsingular  arcs  can  be  very  complex.  For  example,  they  have  shown 
that  for  the  rocket  problem  with  a  bound  on  the  thrust  magnitude,  an  optimal  junction 
between  a  singular  and  an  arc  of  maximum  thrust  must  involve  a  well-defined  infinite 
sequence  of  switches  between  maximum  thrust  and  zero  thrust.  This  infinite  sequence 
of  switches  takes  place  in  a  finite  time.  However,  if  impulses  arc  allowed,  it  is  appar¬ 
ently  possible  to  have  a  simple  junction  between  a  singular  arc  and  an  impulsive  thrust. 
Not  much  work  has  been  done  on  this  problem  and  it  remains  as  a  topic  for  future 
investigation. 

10.3  Singular  Arcs  for  an  Inverse  Square  Field, 

One  of  Lawden's  many  contributions  to  space  trajectory  optimization 
was  his  treatment  of  the  coplanar  singular  arcs,  and  his  analytical  integration  of  the 


time-open  arc  which  has  become  known  as  the  Lawden  spiral.  The  generalized  Legendre- 
Clcbsh  condition  was  subsequently  applied  to  these  arcs  by  Robbins  and  by  Kopp  and  Moyer 
to  yield  the  result  that  the  optimal  thrust  direction  must  point  in  toward  the  center  of  at¬ 
traction  and  that  the  square  of  the  sign  of  the  angle  with  the  horizontal  must  be  less  than 

or  equal  to  1  /ft .  This  result  rules  out  the  Lawden  spiral  as  on  this  singular  arc  the 

/ 

thrust  always  has  a  radial  component  directed  away  from  the  center  of  attraction.  This 
is  a  very  important  result  as  it  indicates  that  for  the  time-open  coplanar  problem  the 
minimizing  solution  can  never  have  a  singular  are.  However,  Robbins  has  shown  that 
for  the  time-fixed  ease  a  significant  prolion  of  the  phase  space  is  filled  by  singular  arcs 
which  do  satisfy  the  generalized  Lcgendre-Clebsli  condition.  Among  those  singular  arcs 


is  one  considered  by  Fraeijs  do  Veubeke  which  corresponds  to  the  fixed-time  angle  open 
case.  This  particular  arc  can  be  integrated  analytically  as  can  the  Lawden  spiral  and  a 
portion  of  this  arc  is  minimizing.  Robbins  uses  an  argument  similar  to  an  argument  in 
the  classical  calculus  of  variations  to  show  that  sufficiently  short,  segments  of  singular 
ai'cs  satisfying  the  generalized  Legendre-Clebsh  condition  are  locally  minimizing.  With 
the  development  of  the  new  sufficiency  conditions  for  singular  arcs,  this  argument  is  no 
longer  necessary  and  it  remains  as  a  topic  for  future  investigation  to  apply  the  additional 
necessary  and  sufficiency  conditions  to  the  singular  arcs  of  an  inverse  square  field.  An¬ 
other  problem  for  future  investigation  is  the  determination  of  composite  extremals  in  an 
inverse  squai-c  field  involving  both  singular  and  nonsingular  arcs. 

The  case  of  singular  arcs  in  three  dimensions  lias  been  treated  by 
Christian  Mnrchal.  Two  of  his  co-workers,  Contonsou  and  Marcc,  have  also  consid¬ 
ered  an  academic  type  of  singular  arc  which  may  occur  in  the  time-open  case.  This 
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peculiar  type  of  extremal  involves  applying  a  finite  acceleration  for  an  iniinitisimal  time 
and  then  coasting  through  an  angle  of  SCO0  before  applying  another  iniinitisimal  pulse. 
These  arcs  do  not  appear  to  be  of  much  interest  because  of  their  somewhat  esoteric 
nature  and  partially  because  they  appear  to  be  nonminimizing. 
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CHAPTER  11:  INTERPLANETARY  TRAJECTORIES 

11. 1  Impulsive  Trajectories 

The  problem  of  calculating  optimum  trajectories  from  one  planet  to  another 
is  an  n-bocly  problem.  For  preliminary  analysis  this  n-body  problem  can  be 
simplified  by  considering  only  the  attraction  of  the  sun  and  of  those  celestial  bodies 
to  which  the  space  craft  makes  a  close  approach  including  the  launch  and  the  arri¬ 
val  planets..  Lawden  has  made  an  approximate  analysis  of  this  n-body  problem  by 
utilizing  what  is  essentially  a  low  order  matched  asymptotic  expansion.  This  type 
of  analysis  has  been  widely  used  for  mission  planning  purposes  and  has  been  found 
to  be  sufficiently  accurate  for  preliminary  design  purposes.  In  this  analysis  the 
radii  of  the  spheres  of  influence  of  the  planets  are  neglected.  The  trajectory  from 
one  planet  to  another  is  calculated  as  a  heliocentric  two  body  trajectory  from  the 
position  of  the  first  planet  to  the  position  of  the  second  planet.  The  relative  veloci¬ 
ties  of  arrival  and  departure  arc  then  regarded  as  the  asymptotic  velocities  on  two 
body  approach  and  departure  hyperbolas  relative  to  the  planets.  The  time  from  one 
planet  to  another  is  based  only  upon  the  heliocentric  two  body  transfe  .  The  time 
inside  the  planetary  spheres  of  influence  is  neglected. 

A  higher  order  analysis  by  Perko  and  Breakwell  of  the  matched  asymptotic 
expansion  shows  that  the  errors  in  calculating  propulsion  requirements  from  this 
model  are  acceptably  small. 

In  calculating  the  primer  vector  for  this  approximate  n-body  trajectory  a 
change  must  be  made  in  the  magnitude  of  the  primer  vector  necessary  to  trigger  an 
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impulse  at  llie  approach  departure  planets.  The  primer  vector  on  the  heliocentric 
trajectory  is  calculated  in  the  normal  fashion  except  that  its  boundary  conditions 
are  somewhat  diflorcnt.  Lawden  has  shown  that  the  magnitude  of  the  primer  vector 
at  the  start  of  the  heliocentric  arc  should  be  smaller  than  unity  by  the  ratio  of  the 
asymptotic  velocity  on  the  initial  departure  hyperbola  to  the  periapsis  velocity  at 
the  beginning  of  this  departure  hyperbola.  A  similar  condition  holds  for  the  arrival 
hyperbola.  These  two  conditions  define  the  magnitude  of  the  primer  vector  at  the 
two  ends  of  the  heliocentric  trajectory.  The  direction  of  the  primer  vector  is  the 
same  ns  the  direction  of  the  relative  velocity  vector  to  the  planets.  The  primer 
vector  on  the  escape  hyperbola  has  the  solution  corresponding  to  the  Lagrange  multi¬ 
plier  for  changing  the  semi-major  axis  or  the  energy  of  the  hyperbola.  The  primer 
vector  is  coincident  with  the  direction  of  the  velocity  vector  and  its  magnitude  is 
proportional  to  the  magnitude  of  the  velocity  vector. 

The  analysis  of  the  heliocentric  trajectory  is  conducted  exactly  as  the  analysis 
of  the  normal  two-body  orbit  except  for  the  changed  condition  on  the  critical  magnitude 
of  the  primer  vector  at  the  launch  and  arrival  planets.  For  an  intermediate  impulse 
given  in  heliocentric  space  the  magnitude  of  the  primer  vector  must  again  have  a 
stationary  maximum  of  unity. 

Two  impulse  transfers  from  one  planet  to  another  are  readily  calculated  by 
solving  Lambert’s  problem  for  given  launch  and  arrival  dates.  If  contours  of  constant 
total  delta  v  are  plotted  as  a  function  of  the  launch  and  arrival  dates  the  well  known 
”pork-chop”  curves  are  obtained.  These  trajectories  may  readily  be  checked  for 
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local  optimality  by  calculating  the  necessary  conditions  on  the  primer  vector  along 
the  trajectory.  During  any  single  opposition  only  portions  of  the  two  impulse  "pork- 


chop"  curves  will  represent  optimal  transfers.  In  many  eases  it  will  be  necessary 
to  go  to  an  additional  mid  course  impulse  or  to  initial  or  final  coasting  periods  in 
order  to  have  an  optimal  transfer  from  one  planet  to  another.  For  some  planets 


such  as  Mars  or  Venus  the  optimal  transfer  during  any  given  opposition  (the  trans¬ 
fer  requiring  the  lowest  total  delta  v  for  any  launch  and  arrival  dales  during  the 
opposition)  will  be  a  two  impulse  transfer.  For  other  planets  such  as  Mercury  an 
optimum  rendezvous  will  generally  require  three  impulses.  An  approximate  first 
order  theory  of  the  optimum  transfer  to  any  planet  during  any  opposition  has  been 
given  by  Marchal.  Numerical  calculations  indicate  that  this  first  order  theory  is 
sufficient  to  be  used  to  determine  initial  conditions  for  numerical  processes  to  find 
the  true  optimum  but  is  not  sufficiently  accurate  to  be  a  good  approximation  to  the 
true  non-linear  minimum.  Numerical  calculations  of  optimum  multiple  impulse 
transfer  from  the  earth  to  other  planets  have  been  carried  out  by  Lion,  Doll  and 
Gobetz.  These  calculations  show  that  mid  course  burns  can  be  used  to  open  up 
launch  windows,  and  for  planets  such  as  Mercury  to  significantly  reduce  the  delta  v 
required  for  rendezvous  missions. 


11.2  Power  Limited  Trajectories 

Many  calculations  have  been  made  of  optimum  low  thrust  trajectories  between 
various  planets.  These  optimum  trajectories  must  at  the  present  time  be  found  by 
numerical  methods  as  none  of  the  approximate  analytic  theories  have  proven  suffi¬ 
ciently  accurate.  Several  fairly  efficient  numerical  techniques  have  been  developed 
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unci  such  trajectories  can  be  routinely  calculated  for  arbitrary  launch  and  arrival 
dates.  'J  here  trajectories  have  been  calculated  l'or  various  assumptions  including 
constant  power,  constant  thrust,  constant  acceleration,  constant  specific  impulse 
and  either  solar  electric  or  nuclear  electric  power  supplies.  The  typical  practice 
is  to  reduce  the  n-body  problem  to  an  approximate  sequence  of  two  body  problems 
as  is  done  for  the  impulsive  case.  An  approximate  theory  of  the  influence  of  the 
departure  planet’s  gravity  on  the  heliocentric  departure  trajectory  has  been 
developed  by  Melbourne  and  by  Edclbnum  by  using  a  low  order  matched  asymptotic 
expansion.  This  analysis  shows  that  a  correction  to  the  heliocentric  trajectory  due 
to  the  planeoeentric  gravity  field  is  necessary,  particularly  for  the  larger  planets. 
The  treatment  of  the  planeoeentric  phases  depends  on  whether  the  energy  to  escape 
from  the  planet  is  provided  by  a  high  thrust  or  a  low  thrust  propulsion  system.  For 
the  case  where  the  space  craft  spirals  away  from  the  planet  under  low  thrust  propul¬ 
sion,  the  time  spent  during  this  escape  maneuver  must  be  considered.  A  refined 
analysis  of  this  problem  has  been  developed  by  Breakwcll  and  Rauch.  For  the  early 
electric  propulsion  systems  it  is  generally  more  desirable  to  use  high  thrust  systems 
to  provide  the  planetary  escape  and  capture  propulsion.  In  these  cases  the  afore¬ 
mentioned  analysis  of  Melbourne  and  Edelbaum  should  be  utilized.  If  curves  of 
constant  fuel  consumption  arc  plotted  against  launch  and  arrival  dates  some  similari¬ 
ties  arc  obtained  with  the  corresponding  results  for  high  thrust  propulsion.  One 
difference  is  that  in  a  low  thrust  case  the  fuel  consumption  generally  decreases 
monotonicnlly  with  the  total  transfer  time.  However,  for  any  given  transfer  time, 
there  is  an  optimum  launch  date  and  the  optimum  launch  dates  for  low  thrust 
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propulsion  lend  to  be  close  to  those  for  high  thrust  propulsion  with  about  the  same 
transfer  times. 

11.3  Stvingbys 

In  one  of  Lawden's  early  papers  he  considered  was  referred  to  as  a  perturba- 
tionnl  maneuver,  where  a  space  craft  on  its  way  to  Mars  would  swing  by  the  moon  in 
order  to  pick  up  some  additional  delta  v,  for  free,  from  the  Moon's  gravitational 
field.  He  found  that  the  delta  v  that  could  be  added  by  the  moon  was  not  vary  large 
and  this  idea  was  neglected  for  about  10  years.  Even  then,  the  idea  was  not  a  new 
one,  as  it  had  long  been  known  that  Jupiter  could  change  the  energy  of  comets  and 
cause  them  to  escape  permanently  from  the  solar  system.  In  the  early  Go's  it  was 
discovered  that  for  round  trip  missions  to  the  near  planets  it  was  often  desirable 
to  use  swingbys  of  Venus  en  route  to  Mars  or  vice  versa.  11  was  also  discovered  that 
the  enormous  gravitational  field  of  Jupiter  could  be  advantageously  used  to  perform 
various  missions  throughout  the  solar  system  including  the  grand  tour  missions, 
close  solar  probes,  and  probes  far  out  of  the  ecliptic  plane.  The  analysis  of  opti¬ 
mum  swingby  missions  falls  within  tic  general  theory  of  optimal  rocket  trajectories 
treated  in  this  monograph.  The  analysis  is  conducted  by  using  matched  asymptotic 
expansions  as  was  the  case  with  planetary  escape  and  arrival  trajectories.  The 
analysis  of  optimal  swingby  missions  can  be  carried  out  both  for  high  thrust  systems 
and  low  thrust  power  limited  rocket  systems.  The  effect  of  the  swingby  is  to  cause 
a  discontinuous  change  in  (he  direction  of  the  primer  vector  at  the  swingby  time  at 
the  two  ends  of  the  heliocentric  approach  and  departure  trajectories  from  the  planet. 
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The  angle  that  the  primer  vector  rotates  through  is  the  same  as  the  angle  through 
which  the  hyperbola  relative  to  the  planet  rotates  its  asymptotic  velocity  vectors. 

It  is  possible  to  treat  both  powered  and  unpowered  swingby  missions  in  this  fashion. 
During  a  powered  swingby  mission  an  impulse  is  applied  at  a  point  close  to  the 
planet.  For  the  powered  swingby  missions  with  high  thrust  it  is  possible  to  use 
the  analysis  of  transfer  between  hyperbolic  asymptotes  treated  in  the  chapter  on 
inverse  square  force  fields.  The  analysis  of  low  thrust  trajectories  during 
swingbys  is  very  similar  to  the  analysis  of  low  thrust  trajectories  during  the 
approach  and  departure  phases  from  the  terminal  planets. 
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CHAPTER  12:  COMBINATION  PROPULSION 


12. 1  Field  Free  Space 

Early  work  on  low  thrust  propulsion  systems  assumed  that  the  low  thrust 
system  would  be  placed  into  orbit  by  a  high  thrust  system  and  that  all  subsequent 
propulsion  would  be  provided  by  the  low  thrust  system.  Later  work  on  the  utili¬ 
zation  of  low  thrust  propulsion  systems  indicated  that  it  was  often  desirable  to 
utilize  high  thrust  in  conjunction  with  low  thrust  for  additional  mission  phases.  It 
lias  been  found  that  combinations  of  high  thrust  and  low  thrust,  propulsion  systems 
can  often  provide  better  performance  than  either  system  when  used  alone.  The 
simplest  such  missions  to  analyse  arc  missions  in  field  free  space.  The  simplest 
of  the  field  free  space  missions  is  the  problem  of  changing  only  velocity.  For  this 
mission  the  payload  with  pure  high  thrust  is  independent  of  the  mission  duration 
while  the  payload  with  pure  low  thrust  increases  monotonioally  with  mission  dura¬ 
tion.  For  mission  times  in  the  vicinity  of  the  time  for  which  the  low  thrust  and  the 
high  thrust  system  have  the  same  payload,  the  combination  of  both  propulsion  sys¬ 
tems  will  provide  higher  payloads. 

The  optimum  combination  load  may  bo  described  as  follows:  a  high  thrust 
impulse  is  first  used  to  accelerate  the  vehicle  to  some  fraction  of  the  final  deshed 
velocity.  Then  the  low  thrust  engine  is  turned  on  and  operates  until  a  later  time 
when  the  power  supply  is  dropped.  After  the  power  supply  is  dropped,  a  second 
high  thrust  impulse  accelerates  only  the  payload  te  the  final  velocity.  For  this  mode 
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of  operation,  the  payload  is  independent  of  the  split  between  the  initial  and  final  im¬ 
pulses.  It  is  possible  to  provide  all  of  the  high  thrust  impulse  initially  and  then  con¬ 
tinue  with  low  thrust  so  that  the  power  supply  is  retained  at  the  end  of  the  mission. 
Alternately,  the  low  thrust  may  be  used  first  after  which  the  power  supply  is  dropped 
and  all  the  high  thrust  is  used  at  the  end.  This  is  a  special  result  which  only  applies 
to  this  particular  mission.  For  this  case  where  the  primer  vector  is  constant  in  mag¬ 
nitude  the  benefits  of  combination  propulsion  are  only  significant  when  the  payload  is 
quite  small. 

A  second  interesting  mission  in  field  free  space  is  a  transfer  between  two  po¬ 
sitions  of  rest  in  field  free  space.  Such  a  transfer  may  be  considered  as  an  approxi¬ 
mation  to  a  fast  transfer  to  the  outer  planets  with  a  nuclear  electric  propulsion  system. 
The  velocities  required  for  fast  transfers  in  this  case  arc  so  high  that  the  influence  of 
gravitational  fields  and  terminal  velocities  may  be  neglected.  The  optimum  mode 
with  pure  high  thrust  is  to  provide  an  initial  impulse,  coast,  for  a  given  time  and  then 

decelerate  at  the  target  position  back  to  zero  velocity  with  a  second  impulse.  In  this 
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case  the  high  thrust  payload  is  a  function  of  the  flight  time  and  increases  monotonically 
with  the  flight  time  as  does  the  low  thrust  payload.  For  short  transfer  times  the  pay- 
load  with  the  high  thrust  system  will  be  larger  than  the  payload  of  the  low  thrust  sys¬ 
tems  while  for  long  transfer  times  the  payload  with  the  low  thrust  system  will  be  larg¬ 
er.  Once  again  combination  propulsion  is  advantageous  in  the  region  where  both  sys¬ 
tems  have  about  the  same  performance.  The  combination  mode  for  this  mission  is  to 
first  provide  an  initial  velocity  impulse  ,  and  then  turn  on  the  low  thrust  which  has  an 
acceleration  magnitude  that  decreases  linearly  with  time.  The  acceleration  magnitude 
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passer,  through  zero  and  then  increases  to  its  initial  acceleration  magnitude.  A  sec¬ 
ond  and  equal  high  thrust  impulse  is  then  provided  after  the  power  supply  has  been 
dropped.  Once  again  the  combination  propulsion  system  can  provide  appreciable 
payload  increases  where  the  payload  of  either  system  alone  is  quite  small  and  about 
equal  to  the  payload  of  the  other  system. 

12.2  inverse  Squ are  Force  Fields 

In  more  complicated  gravitational  fields  such  as  inverse  square  fields  the 
difference  between  the  maximum  and  the  average  values  of  the  primer  vector  along 
the  optimum  trajectory  maybe  much  greater  than  for  the  eases  in  field  free  space. 

By  utilizing  high  thrust  systems  in  those  regions  where  the  primer  vector  has  its 
maximum  values,  much  greater  bene  Cits  from  combination  operation  arc  possible 
than  for  the  missions  in  field  free  space.  If  the  inverse  square  field  is  a  strong  field 
then  the  analysis  of  Chapter  7  may  be  used  for  the  low  thrust  phases  of  the  mission. 

The  high  thrust  phases  may  use  the  analysis  in  Chapter  9.  Since  the  high  thrust  will 
be  used  impulsively  the  low  thrust  system  may  be  assumed  to  be  on  all  the  time  except 
possibly  at  the  final  impulse  before  which  the  power  supply  will  be  dropped.  If  there 
is  more  than  one  high  thrust  impulse  it  might  be  advantageous  to  drop  portions  of  the 
low  thrust  power  supply  during  the  mission.  This  will  also  be  true  if  the  payload  for 
the  low  thrust  portion  of  the  system  becomes  very  small  so  that  optimal  staging  of 
the  low  thrust  power  supply  becomes  desirable.  The  optimal  mission  modes  in  inverse 
square  fields  become  very  complicated  and  in  different  regions  different  combinations 
and  sequences  of  low  thrust  and  high  thrust  become  optimal.  It  is  necessary  to  assume 
different  modes  and  check  the  local  optimality  by  calculating  the  primer  vector  in  order 


to  determine  the  local  optima.  These  local  optima  must  then  be  compared  by  deter¬ 
mining  which  of  the  different  modes  tends  to  lead  to  the  absolute  minimum.  The  mis¬ 
sion  that  has  been  analysed  in  most  detail  is  escape  from  an  inverse  square  field  with 
a  given  hyperbolic  velocity.  The  optimum  mode  tor  this  mission  has  generally  been 
found  to  be  an  initial  period  of  low  thrust  which  increases  the  semi-major  axis  of  the 
initial  circular  orbit  and  also  increases  the  eccentricity.  At  the  end  of  the  first  low 
thrust  phase  the  vehicle  will  be  in  an  elliptic  orbit,  usually  with  a  minimal  allowable 
periapsis  radius  and  a  faii'ly  high  apoopsis  radius.  At  this  point  a  high  thrust  impulse 
is  applied  at  periapsis  to  accelerate  the  vehicle  beyond  escape  energy.  A  second  pe¬ 
riod  of  low  thrust  may  then  be  used  if  the  final  desired  hyperbolic  velocity  is  fairly 
large. 
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For  transfer  between  co-planer,  co-axial  elipscs  the  optimum  mode  may  re¬ 
quire  from  zero  to  two  initial  impulses.  A  second  impulse,  if  it  occurs,  being  given 
,  at  the  second  apsis  of  the  first  transfer  clipse.  And  there  may  also  be  from  zero  to 
two  terminal  impulses.  Usually,  there  will  be  no  more  than  two  high  thrust  impulses 
in  combination  with  the  low  thrust  system,  although,  occasionally  three  impulses  may 
be  required.  For  the  cases  with  two  high  thrust  impulses,  both  impulses  might  occur 
at  the  beginning,  one  may  occur  at  the  beginning  and  one  at  the  end,  or  both  impulses 
may  occur  at  the  end  depending  upon  the  particular  initial  and  final  orbits. 
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PROBLEMS 


1.  A  body  moves  in  a  central  force  field  where  the  force  is  directed  towards 
the  origin  and  is  proportional  to  the  distance  from  the  origin.  What  shape 
does  an  unpowered  orbit  have  in  such  a  force  field?  Mow  many  periapses 
and  apoapsides  tire  there  in  one  revolution?  What  is  die  period  of  the  orbit? 
(Rectangular  coordinates  are  suggested  for  this  and  the  following  two  prob- 
lcms.) 

2.  Assume  that  (he  optimum  impulsive  transfer  between  a  circular  orbit  and 
a  coplanar  elliptic  orbit  in  the  above  force  field  requires  two  tangential 
impulses  at.  the  apsides.  What  is  the  minimum  A  V  required  for  transfer 
from  a  circular  (o  an  elliptic  orbit?  Does  the  transfer  to  the  apoapsis  or 
to  the  peri  apsis  require  less  fuel?  Both  the  periapsis  and  the  apoapsis  of 
the  ellipse  may  be  either  larger  or  smaller  than  the  radius  of  the  circular 
orbit. 

3.  Derive  the  optimum  thrust  program  for  a  power-Jimited  vehicle  in  the 
linear  central  force  field  of  Problem  1.  Express  this  program  in  terms 

of  the  initial  values  of  the  Lagrange  multipliers  and  the  time  by  integrating 
the  Euler-Lagrange  equations. 

4.  A  power-limited  vehicle  is  moving  along  the  x  axis  towards  the  y  axis  in 
field-free  space.  What  is  the  minimum  J  to  transfer  from  velocity  u°  and 
position  x°  to  a  final  velocity  directed  along  the  y  axis  in  the  time  i  \  ? 
What  is  the  optimum  value  of  y  at  U  for  this  J?  What  is  the  optimum  value 
of  x°? 

5.  Synethesize  the  optimal  control  for  transfer  from  any  state  to  the  final 
state  of  problem  4. 
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